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ABSTRACT 

^  The  minliaum  information  approach  to  active  control  of  struc¬ 
tural  systems  seeks  inherently  robust  designs  by  use  of  mean- 
square  optimization  conjoined  with  a  stochastic  system  model 
which  presvones  as  little  as  possible  regarding  a  priori  informa¬ 
tion  on  modal  parameter  statistics.  This  report  extends  earlier 
results  for  the  regulator  problem  to  the  case  of  full-order 
dynconic  compensation  with  nonsingular  observation  noise.  Opti¬ 
mality  conditions  along  with  sufficient  conditions  for  existence 
and  uniqueness  of  solutions  and  for  closed- loop  stochastic 
stcd^ility  are  presented.  Results  concerning  asymptotic  proper¬ 
ties  for  large  uncertainty  levels  are  also  given.  Numerical 
results  for  various  simple  examples  indicate  improved  robustness 
properties  over  standard  LQG  designs  and  suggest  the  possibility 
that,  under  the  minimum  information  stochastic  approach,  the 
burden  of  design  ccanputation  may  be  reduced  to  that  associated 
with  the  relatively  well  known  or^coherenf^odes. 
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MEAN-SQUARE  OPTIMAL,  FULL-ORDER  COMPENSATION  OF 
STRUCTURAL  SYSTEMS  WITH  UNCERTAIN  PARAMETERS 


1.  INTRODUCTION 

Many  techniques  advanced  for  the  active  control  of  highly 
flexible  structural  systems  implicitly  assume  the  existence  of 
a  large  order  "verification"  model  in  which  the  values  of  struc¬ 
tural  parameters  are  precisely  known.  Where  the  critically 
important  issue  of  parameter  uncertainties  (arising  from  inher¬ 
ent  limitations  in  the  analysis  of  high  order  structural  modes) 
is  specifically  addressed  within  a  computationally  tractable 
approach,  the  formulation  entails  either  (a)  an  essentially 

geometric  (in  the  linear  algebraic  sense)  approach  to  sensitiv- 
* 

ity  reduction,  (b)  ad  hoc  combination  of  inherently  robust  con- 

** 

trols  with  more  standard  LQ  design  techniques  and/or  (c)  the 
a  posteriori  verification  of  robustness  by  use  of  methods 
entirely  extrinsic  to  the  design  process.^ 

In  contrast,  it  may  be  argued  that  although  verification 
models  are  precise,  they  are  necessarily  false;  that  the  struc¬ 
ture  must  be  regarded  as  a  stochastically  parcunetered  mechani¬ 
cal  system  (for  which  limited  a  priori  statistical  data  is 
unavailable) ;  and  that  design  optimization  must  proceed  from  a 
measure  of  performance  defined  over  the  entire  parameter 

*  For  example,  the  work  of  Sesak  [1,2]  as  reinterpreted  by 
Coradetti  [3]  actually  involves  static  decoupling  of  sensitivity 
derivatives  and/or  modal  coordinates  for  model  order  reduction. 

**  e.g. ,  the  rate-output  feedback  of  Balas  [4]  or  low-authority 
control  of  Aubrun  [5]  would  be  charged  with  the  control  of  high- 
order,  poorly- known  modes  while  providing  enhanced  stabilization 
of  a  "high  authority"  LQ  control  designed  for  a  reduced  order 
model  of  the  low  frequency  modes. 

t  Using,  for  example,  general  multivariable  measures  of  robust¬ 
ness  along  the  lines  of  Ref.  [6]. 
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statistical  ensemble.  Very  recently,  such  a  stochastic  design 
approach  was  outlined  and  specifically  elaborated  for  the  special 
case  of  full  state  feedback  regulation  of  structural  systems  with 
a  priori  uncertainties  in  the  modal  frequencies  [7,8]. 

Entailing,  as  it  does,  a  fundcunental  revision  of  the  tradi¬ 
tional,  deterministic  approach  to  dynamical  modelling,  the  ear¬ 
lier  developments  [7-9]  must  be  reviewed  here.  To  retain  the 
spirit  of  linear-quadratic  optimization,  we  choose  as  the  per- 
form£mce  measure  the  average  of  a  quadratic  functional  of  state 
and  controls  over  the  parameter  ensemble.  Secondly,  as  a  com¬ 
plete  ^pirical  specification  of  parameter  statistics  is  never 
provided  in  practice,  the  design  approach  must  accept  the  kind 
of  severely  limited  statistical  data  that  is  actually  available. 

To  avoid  ad  hoc  assumptions,  the  full  probability  assignment 
required  for  determination  of  the  mean-square  optimal  control 
must  be  consistent  with  the  data  on  hand  but  maximally  unpre- 
sumptive  with  regard  to  unavailable  data.  This  is  accomplished 
by  resort  to  a  maximum  entropy  principle  (Jaynes'  principle). 

Furthermore,  to  achieve  particular  simplicity  and  design 
conservatism,  we  acknowledge  as  "available"  the  minimum  possible 
a  priori  statistical  data  needed  to  Induce  a  complete  probabil¬ 
ity  assignment  while  preserving  fidelity  of  the  overall  model  at 
high  levels  of  uncertainty  (or  for  high  order  modes) . 

Thus,  the  approach  is  "minimum  information"  in  two  respects — 
first,  -  because  we  acknowledge  as  available  the  bare  minimum  of 
a  priori  data  and,  secondly,  because  we  induce  a  full  probabil¬ 
ity  model  from  this  acknowledged  data  by  use  of  a  minimum  infor¬ 
mation  (i.e. ,  maximum  entropy)  principle. 
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With  uncertainties  only  in  the  open  loop  frequencies,  we 
identified  the  minimum  data  set  by  examining  the  phenomenology 
of  frequency  uncertainties  as  reflected  in  the  mean  (the  par¬ 
ameter  ensemble  averaged)  response,  the  covariance  matrix  and 
the  "expected  cost"  matrix  (i.e. ,  the  covariance  of  the  co-state). 
The  principal  effects  of  modal  frequency  uncertainty  are  the 
introduction  of  a  spurious  deunping  (the  "decorrelation  damping") 
into  the  mean  response  and  the  suppression  of  cross-correlation 
auaong  distinct  modes.  Essential  to  the  proper  modelling  of  such 
qualitative  features  are  the  decorrelation  deunping  time  con¬ 
stants,  termed  the  "modal  decorrelation  times".  These  are  gen¬ 
erally  inversely  proportional  to  modal  frequency  standard  devia¬ 
tions  and  constitute  fundamental,  albeit  unconventional,  measures 
of  frequency  uncertainty. 

Acknowledgement  of  only  the  mean  (or  nominal)  values  of 
modal  frequencies  and  the  modal  decorrelation  times  as  available 
data  induces  a  white  parameter  statistical  model  which  reduces 
the  optimization  problen  to  solution  of  a  modified  Riccati 
equation  (the  "stochastic  Riccati  equation")  for  the  expected 
cost  matrix.  Under  mild  restrictions,  this  possesses  a  unique 
positive  semi-definite  solution  which  guarantees  closed-loop 
stochastic  stability.  Thus  the  proposed  approach  reduces  the 
need  for  design  iteration  to  achieve  robust  stability. 

The  most  significant  aspect  of  the  stochastic  Riccati 
equation  is  the  character  of  its  steady-state  solutions  for 
large  uncertainties  [9].  If  the  uncertainties  in  all  open-loop 
frequencies  increase  without  bound  (i.e.,  all  decorrelation 
times  approach  zero)  the  expected  cost  reduces  to  a  diagonal 
matrix  whose  elements  are  independent  of  modal  frequency  statis¬ 
tics  and  are  given  by  simple  analytical  expressions.  This 
asymptotic  solution  gives  rise  to  a  velocity  feedback  control  law 
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which  is  stable  regardless  of  the  values  of  modal  frequencies  or 
damping  ratios.  Thus  a  greater  degree  of  robustness  is  obtained 
than  was  originally  sought. 

In  the  more  typical  case  in  which  the  frequencies  of  low 
order  structural  modes  are  relatively  well  known  while  modelling 
accuracy  deteriorates  for  the  high-order  modes,  the  stochastic 
Riccati  equation  automatically  produces  a  velocity  feedback  con¬ 
trol  (of  the  asymptotic  form)  for  the  high-order,  poorly-known 
modes.  At  the  same  time,  for  low-order  modes  having  small  uncer¬ 
tainties,  the  control  closely  resembles  the  deterministic  plant 
solution.  In  other  words,  a  "high  authority",  essentially  deter¬ 
ministic  control  for  well-known  modes  and  a  "low  authority" 
velocity  feedback  control  for  relatively  uncertain  modes  natu¬ 
rally  anerge  as  limiting  regimes  of  a  jlobal  control  law  which 
is  guaranteed  to  be  stable  over  the  parameter  ensemble. 

This  general  behavior  has  immediate  consequences  for  the 
computational  effort  required  for  high  order  systems.  In  brief, 
the  computational  burden  required  for  determination  of  the  con¬ 
trol  gain  is  mainly  associated  with  the  relatively  few  weJ.l-known 
("coherent")  modes.  Provided  that  the  dimension  of  the  coherent 
system  is  moderate,  the  stochastic  Riccati  equation  is  amenable 
to  numerical  solution  of  acceptable  accuracy  for  systems  of 
arbitrary  order. 

The  above  qualitative  features  provide  strong  motivation  for 
extension  of  the  formulation  to  linear,  dynamic  compensation. 

Here  we  consider  a  structural  system  as  described  by  a  finite 
number  of  its  "all-elastic"  normal  modes  with  uncertainties  in 
the  open- loop  frequencies.  We  remove,  however,  the  restriction 
of  previous  work  to  full-state  feedback  and  consider  dynamic  out¬ 
put  feedback  ccmpensation.  Chapter  2  sets  forth  the  minimum 
Information  stochastic  modelling  approach  and  derives  the 
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stochastic  Lyapunov  equation  which  must  be  appended  to  the  mean- 
square  optimization  problem  as  a  constraint.  Various  properties 
of  this  equation  are  related  to  closed-loop  stochastic  stability 
and  the  existence  of  steady  state,  constant  gain  controls.  In 
Chapter  3  we  specialize  to  full-order  dynamic  compensation,  and 
derive  explicit  stationary  conditions.  Existence  and  uniqueness 
of  solutions  to  these  conditions  as  well  as  their  asymptotic 
properties  for  large  levels  of  modal  frequency  uncertainty  are 
explored.  Chapter  4  concludes  by  a  presentation  of  computational 
procedures  and  various  numerical  examples.  From  these  theoreti¬ 
cal  developments  and  supporting  numerical  results,  it  will  be 
seen  that  the  desirable  features  of  the  stochastic  design  approach 
as  applied  to  full-state  feedback  regulation  are  indeed  retained 
in  this  less  idealized  setting. 
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THE  MINIMUM  INFORMATION  FORMULATION  OF  THE  OPTIMIZATION 
PROBLEM 


2.1  Problem  Statement 

As  in  the  previous  work  [7],  we  consider  the  control  of  a 
linear  elastic  structure  subject  to  small  deformations  and  no 
rigid  body  degrees  of  freedom.  Modifications  needed  to  include 
rigid  body  modes  are  straightforward  and  need  not  be  treated  in 
this  initial  development.  Retaining  n  normal  mode  coordinates 
in  the  system  model,  the  state-space  form  of  the  equations  of 
motion  may  be  written: 


X  =  (A  +  a(t))x  +  Buj^  + 


n  =  Cx  +  w. 


■r  _2nx2n  _  _2nx£,  p 

AeR  ,BeR  ,u,  eR,neR 


where  x  is  the  vector  of  modal  coordinates  and  velocities  with 
its  odd  indexed  elements  representing  modal  displacements  and 
the  adjacent  even  indexed  elements  giving  the  corresponding 
modal  velocities,  w^^  is  a  white  disturbance  noise  with  inten¬ 
sity  2  0,  Uj^  the  control  input  and  B  the  input  map: 


21  "22 


41  ^42 


where  the  non- zero  elements  are  proportional  to  the  normal  mode 
shapes  at  I  actuator  locations.  n  is  the  vector  of  sensor  out¬ 
puts  with  output  map  C  and  observation  noise,  W2«  We  assume 
that  W2  is  independent  of  and  has  a  nonsingular  intensity 


matrix  V2. 


X  is  the  nominal  or  mean  value  of  the  system  map: 


A  =  block-diag 
ksl, . . . ,n 


where  the  k=l,...,n  are  the  nominal  design  values  of  the 
modal  frequencies  and  the  are  the  modal  damping  ratios 
representing  inherent  structural  damping.  It  is  assumed 
throughout  that: 


Finally,  a(t)  is  the  random  portion  (assumed  zero  mean)  of  the 
system  map  representing  possible  statistical  variation  of  the 
modal  frequencies. 

As  a  preliminary  step,  it  is  convenient  to  express  the 
above  relations  in  the  eigen-basis  of  A.  In  view  of  the  assump 
tion  of  small  damping  we  may  simplify  this  process  by  introduc¬ 
ing  the  resonant  approximation  for  A: 

_  ^ 

A  =  block-diag  ^2  —  ^ 

k“l  ,...,n  _—  —  n  . 

so-called  because  the  difference  between  damped  and  undamped 
natural  frequencies  is  neglected.  With  this  replacement,  the 
eigenvector  matrix  of  A  is: 

r  1  1 

♦  =  block-diag  _  _  ( 

k“l , . . . ,  n  _  — m)j^ _ 


Then,  defining: 


=  diag  { ,Wj^(-i-n2^)  / 


a. 


v(t)  ^  a  ♦ 

B  =  B 

~  A  *-l 

i 

^  *-l  T7  *-lH 

-9  Vj^  $ 

A  o  * 

Y  =  C  $ 


the  equations  of  motion  written  in  terms  of  the  state- space 
modal  coordinate  vector: 


C  ^  X 


assume  the  form: 


C  =  uC  +  6  +  Wj^ 


n  =  yC  +  w. 


y  =  y  +  v(t) 


Furthermore,  we  suppose  that  a  priori  uncertainty  exists 
only  in  the  open-loop  frequencies  so  that: 


v(t) 


(10) 


«  diag  {iiro  (y.)  6.  (t)  } 

k-l,...,2n  ^  ^ 

where  Iin(...)  denotes  the  Imaginary  part.  The  6j^(t)  ,  k=l,...,2n, 
are  assvimed  real  valued,  zero-mean,  stationary  random  processes 
in  time  and  statistically  independent  of  w^  and  ^2’ 

Obviously,  a  general  treatment  would  require  that  all  sys¬ 
tem  maps  be  random.  However,  the  above  restriction  offers  an 
appropriately  simple  point  of  departure  and  permits  relatively 
easy  interpretation.  Moreover,  as  will  ultimately  be  seen,  the 
very  special  problem  considered  here  still  exhibits  important 
features. 

With  the  above  restrictions,  suppose  the  control  to  be  pro¬ 
vided  by  interconnection  of  (9)  with  a  fixed  order  dynamic  com¬ 
pensator.  Specifically,  the  controlled  system  equations  take 
the  form: 


e 

C  = 

yC  + 

0  Uj^  +  Wj^ 

;  5 

^  /^2n 
e  C 

a.  V 

!  (11) 

q  = 

aq  + 

"2 

;  q 

e  C 

b.  ) 

where 

^  -icq 

n  »  y5  +  Wj 
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( 


Intensity 


0  1 


;  2  0,  V2  >  0 


and  where  q  is  the  compensator  state  and  k,£  are  the  control 
gains. 

As  a  measure  of  performance  of  system  (11) ,  we  take  the 
quadratic  functional: 


7  =  E 


dt  +  u”  R2Uj^1 

o 

lx  I 


R2  >  0  ,  e  R 


r^^  *  *  Rj^  ♦ 


R^  i  0  I  c  R 


2nx2n 


a. 


b. 


c. 


d. 


C 


where  the  averaging  operation  includes  the  parcuneter  ensemble 
(i.e.,  in  this  case  the  ensemble  of  the  6j^(t)).  Note  that  term¬ 
inal  state  weighting  is  ignored  in  (14)  since  its  inclusion 
would  entail  only  superficial  modification  of  the  following 
results. 

The  problem  is  to  determine  k  and  f  in  (12)  (with  appro¬ 
priate  choice  of  a)  to  minimize  J.  It  is  expeditious,  at  this 
point  to  recast  this  problem  in  terms  of  the  augmented  state: 


in  which  case  (11)  through  (14)  become: 


min:  1  »  E[/  ^  dt  x“  R'X] 

K ,  f ,  o  ^o 

X  =  [A'  +  V*  -F'  ]  X  +  W(t) 


where 


and 


0 

a  . 


(16) 


(17) 


(18) 


a. 


b. 


c.  >.  (19) 


d. 


e. 


Finally,  after  straightforward  manipulation,  we  can 
restate  the  optimization  problem  in  terms  of  the  augmented 
co-state  matrix  as  follows: 


min:  J 

,a 


(20) 


=  tr  / 


dt  P' 


(t)  V 


f  .  E[P'l  e 


a. 


-P'  =  [A*  +  V  -F']“  P'  +  P' [A'  +  V  -F'l  +  R'  b.  >  (21) 


P'  (tj^)  =  0 


c. 


where  the  averaging  in  (21. a)  now  extends  only  over  the  parameter 
ensemble.  The  quantity  P'  (which  we  term  the  expected  cost 
matrix)  is  the  covariance  of  the  augmented  co-state  and  gives  a 
direct  measure  of  mean-square  performance  by  virtue  of  (20) . 

From  the  assumptions  made  concerning  W(t) ,  the  disturbance 
noise  is  a  differential  process  (with  unbounded  variation)  and 
care  is  needed  in  the  precise  specification  of  the  meaning  of 
(17).  Such  a  specification  should  also  allow  the  parameter 
noise,  v' ,  to  be  a  differential  process.  To  answer  this  need, 
we  shall  talce  (17)  to  mean  the  truth  with  probability  one  of  the 
equality: 


/^h(t)  [-dX(t)  +  [A'-F' lX(t)dt  +  dv'(t)X  +  dW(t)  ]  =  0  (22. a) 

a 


for  all  t  ^  a  <  b  ^  t,  and  any  continuous  matrix  function,  h(t)  ; 
o  X 

where  the  above  stochastic  integrals  are  interpreted  according  to 
Stratonovich  [10,11]  and  Wong  and  Zalcai  [12]  : 


(22 


f\(t)  dv’(t)  X(t)  - 

A 


l.i.m.  2h(t.)  [V  (t.)  -V  (t._i)]  [%(x(tj^)  +  X(ti_i))] 
e-^0  i 


where  is  a  partition  of  the  interval  [a,b]  and 


£  ^  max  (t^ 

i 


I^i-l/^il 


(23) 


Interpretation  (22)  is  the  appropriate  one  for  our  applica¬ 
tion  since,  in  writing  (17)  we  have  in  mind  a  mechanical  system 
whose  parameters  are  perturbed  by  a  noise  of  finite  total  power. 
Adopting  the  more  familiar  Ito  differential  for  (17)  would 
ignore  this  fact.  However,  in  the  case  in  which  v*  is  white, 

(22)  yields  results  corresponding  to  a  bandpass  pareuneter  noise 
in  the  limit  as  the  passband  approaches  infinity. 

Thus,  with  (22) ,  v*  may  be  treated  as  a  process  of  bounded 
vauriation  (almost  everywhere)  and  we  may  state  various  formal 
results  for  the  system  response  in  terms  of  the  transition 
matrix  of  (A*  -F'  +  V)  as  follows: 

Theorem  1 

Suppose  that  V (t)  is  a  stationary  zero  mean  random  matrix 
process.  Define  an  increment  in  the  nonstationary  process  W  by: 

-  W'(t.)  -W'(t,)  -  fj  dx  v'(t)  ;  t.,  i  t. 


W  (t,  ,t.,) 


(2 


and  assume  that  W'(t^,t2)  possesses  joint  moments  of  all  orders 
for  all  finite  t2-tj^.  Further,  suppose  that  F'  (t)  is  bounded 
and  continuous.  Then  under  interpretation  (22)  and  in  a  prob¬ 
ability-one  sense: 

A.  The  transition  matrix,  <f'{t,T),  for  system  (17)  is  given  by 


<!>' (t,  t) 

\«here 


2  (t,T)  ;  t  >  T 

k=0  k 


a.  ^ 


<|)Q(t,T)  =  exp  [A*  (t-x)  +  W'(T,t)] 

T  T 

<)>j^(t,T)  =  (-1)^/!^  dXj^  dT2.../^^"^  dTj^ 

(k>0) 

X  Uo(t,Tj^)  F'  (Xj^)  (j»^(Tj^,T2)  F'(T2)... 


b. 

r  (2 


c.  y 


and  where  the  integrals  are  the  usual  Riemann-Stielt jes  sums  and 
extend  over  the  left  semi-closed  Intervals. 

B.  (t>'(t,T)  is  almost  everywhere  continuous  and  its  first  and 
second  moments  are  continuous  and  differentiable  in  both 
arguments. 


C. 


Egs.  (21.b,c)  possess  the  unique,  positive  semi-definite 
solution: 


P'  (t) 


dT  \|<(t,T)  ;  t  e 


(26 


where,  for  t  ^  t,  6s  T-t; 

i|;(t,T)  =  (^•“(t+6,t)  ij)(t+6,T)  (ji'(t+6,t)  a.  j 

!  (27 

ij;(T,T)  =  R'  (t)  b.  ^ 

D.  ip(t,T)  ^  Eli|;(t,T)];  t  e  continuous  and  dif¬ 

ferentiable  in  t  and  t. 

The  above  are  entirely  analogous  to  the  results  of  Theorem  1 
of  [7]  and  the  proof  may  be  omitted  here. 

In  the  case  in  which  v'  is  a  differential  process,  we  have 
an  alternative  formulation.  In  such  a  case,  by  virtue  of  (22) , 
system  (17)  possesses  the  I to  differential: 

d  X'(t)  =  (A'  -F'(t)  +  hi’)  X(t)  dt 

+  dv'  (t)  X(t)  +  dW(t) 

J'  =  lim  i  E[W'^(t,t+A)  ] 

A+0 


where  ^sJ*  is  the  so-called  Stratonovich  correction. 


2. 2  The  Minimum  Information  Model  -  Derivation  of  the 
Stochastic  Lyapunov  Equation 

As  may  be  seen  from  Theorem  1,  a  complete  specification  of 
the  statistical  structure  of  open-loop  frequency  deviations  per¬ 
mits  explicit  determination  of  P' .  However,  a  complete  proba¬ 
bility  model  of  the  parcuneters  based  upon  empirical  determina¬ 
tions  can  never  be  provided  in  practice  and  we  are  faced  with 
limited  available  data  on  parameter  statistics.  To  induce  a 
ccxnplete  probability  model  uniquely  from  the  available  data  we 
define  the  desired  probability  assignment  as  the  one  which,  under 
the  constraints  imposed  by  available  data,  is  maximally  noncom- 
ittal  with  regard  to  unavailable  data,  i.e.,  maximizes  the 
entropy  of  the  underlying  processes. 

More  specifically,  we  may  introduce  measures  of  information 
reposed  in  the  statistics  of  relative  deviations  of  modal  fre¬ 
quencies  in  the  following  way.  First,  note  that  only  increments 
in  the  non- stationary  processes  6j^(0,t);  defined  by 


A  ^2 

=  6j^(0,t2)  -6j^(0,t^)  = 

t2  ^  tj^  ;  k  =  1,  .  .  . ,  2n 


(29) 


actually  enter  into  (26)  and  (27)  and  thus  need  to  be  considered, 
(k) 

Define  {t  ' }  as  an  arbitrary  division  of  the  real  line  with 
m 

m  =  0,1,..., Nj^;  Nj^  finite  and  k  e  where  N  comprises  a  set  of 
distinct  integers  in  the  range  1  to  2n.  For  notational  conven¬ 
ience  let  be  the  totality  of  such  divisions.  Then  a 

measure  of  ^  information  contained  in  the  finite  set  of 
increments : 
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(30) 


k  e  N 


is  given  by  the  relative  entropy  of  these  increments. 


H[6;{t^^^^M]  =  -/  da  P[6;{t^^^^^  }]  In  P  [6;  { }]  (31) 

“k  “k  “k 

where  da  is  the  volume  element  in  the  sample  space  of  (30)  and 

P[5;{tA  }]  is  the  joint  probability  density  of  the  increments 

"k 

(30) .  Although  it  is  possible  to  define  a  measure  of  entropy 
for  the  stochastic  system  as  a  whole,  the  measures  (31)  defined 
directly  on  the  ®j^(t2^rt2)  will  suffice  for  present  purposes. 

Now  in  practice,  we  may  suppose  that  specific  numerical 
values  may  be  assigned  to  various  statistics  of  the  6^(t, ,t,) 
(defined  as  functionals  of  Pl6;{t^  '>])  based  on  empirical 

determinations.  This  constitutes  ^  the  "available  data"  of  the 
problem.  To  avoid  ad  hoc  assumptions  on  the  probability  distri- 
butions  of  frequency  uncertainties,  we  choose  P[6;{t^  Ml  to 
maximize  a  measure  of  our  ignorance  of  the  increments^ (30)  in 
the  light  of  the  available  data.  In  other  words,  given  con¬ 
straints  on  P  [6;  >1  implied  by  the  available  data,  we 

determine  ^  P  [6;  }]  to  maximize  H  [6;  }]  of 

(31)  for  all  choices  of  JV,  n!^  and  {t^^^  ;  m  e  [0,N,  ],  ^  k  e  ff}. 

jc  m  Jc 

The  idea  of  employing  a  statistical  model  which  is  maxi¬ 
mally  unpresumptive  with  regard  to  parameter  data  can  be  carried 
still  further.  We  may  choose  to  acknowledge  as  available  a  data 
set  which  is  essential  to  the  proper  modelling  of  open-loop 
statistical  response  and  is  just  sufficient  to  induce  a  well- 
defined  probability  assignment  via  a  maximum  entropy  principle. 
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The  discussion  of  section  (3)  of  [13]  established  various  quali¬ 
tative  features  of  the  open-loop  system  response  (in  particular 
the  mean,  covariance  and  expected  cost) ,  and  introduced  the  con¬ 
cepts  of  "decorrelation  damping",  "coherence  limit"  and  "incoher 
ent  range".  It  was  concluded  that,  at  the  very  least,  any 
approximating  probability  model  of  the  6j^(0,  t)  ;k=l, . . . ,  2n  should 
preserve  the  time  scales  of  decorrelation  dcunping,  provide  a  cor 
rect  estimate  of  the  coherence  limit  and  satisfy  the  bound  given 
by  (19. b)  of  [13]  for  the  cross-correlations  of  high  order  modes 
This  is  possible  only  if  the  "modal  decorrelation  times",  T  • 


Tj^  ^  (  llm(yj^) 


=  dt|E[exp  iim(yj^)  dx  6j^(t)]1^ 


(3 


k  =  1, . . . , 2n 


'^2m  “  '^2m-l  '  ™ 


are  admitted  as  fundamental  data.  In  essence,  numerical  values 
assigned  to  the  Tj^  establish  the  scales  of  frequency  deviations 
relative  to  the  remaining  time  scales  of  the  problem. 

Thus,  we  propose  to  acknowledge  only  the  Tj^ (k=l, .  . . ,  2n)  as 
the  "available"  data.  It  remains  to  determine  the  probability 

assignment  which  maximizes  the  entropy  (31)  for  all  N,  N.  , 

(k)  ^ 

{t  ;  m  e  [0,N,],  keN}  subject  to  the  constraints  implied  by 
m  jc 

(32)  .  We  term  the  resulting  probability  model  the  "maximum 
entropy  probability  assignment  induced  by  the  data"  (32) . 

Here  we  repeat  the  answer  to  this  problem  given  in  Theorem 
3  of  [7] : 


Theorem  2 


Assuming  that  the  processes  (29)  possess  finite  variances 
for  all  t^,t2  e(0,<»)  and  stationary  increments,  the  maximum 
entropy  probability  assignment  induced  by  (32)  is  the  one  under 
which  the  6j^(0,  t)  ;k=l, . . . ,  2n  are  independent  Wiener-Levy  pro¬ 
cesses  with  intensities  • 


Er6r  (0,t)]  =  I.  |t|/|lmy. 


^m  “  ^2m-l  '  ro  -  wn 


It  must  be  noted  that  if  the  available  data  encompasses 
more  than  specification  of  the  decorrelation  times,  the  result¬ 
ing  maximum  values  of  the  infoinnation  measures  (31)  are  smaller 
than  the  values  corresponding  to  the  above  probability  assign¬ 
ment.  On  the  other  hand,  if  the  available  data  omits  some  of 
the  decorrelation  times,  a  maximum  entropy  probability  a^sign-^ 
ment  for  which  (29)  possess  finite  variances  for  finits:  t  does 
not  exist.  This  arises  because  (31)  monotonically  increases 
with  the  noise  intensities  so  that,  loosely  speaking,  the  maxi¬ 
mum  entropy  model  involves  white  noise  of  unbounded  intensity. 
Thus,  in  this  sense,  the  decorrelation  times  constitute  the  min¬ 
imum  data  required  to  induce  any  "reasonable"  maximum  entropy 
probability  assignment. 

Most  importantly,  the  ^]^(t)  are  modelled  as  white  noise  so 
that  posterior  learning  is  impossible  and  the  stochastic  control 
problem  is  nondual.  In  consequence,  the  white  parameter  uncer¬ 
tainty  model  provides  a  worst  case  situation  from  the  point  of 
view  of  parcuaeter  identification.  Indeed  the  model  may  be  used 


K 


to  determine  performance  degradation  due  to  parameter  uncertainty 
and  to  assess  the  need  for  identification  and  adaptive  algorithms. 

With  this  maximally  unpresumptive  statistical  model  we  are 
in  a  position  to  determine  a  single  closed  equation  for  P'  and 
restate  the  optimization  problem  as  follows; 

Theorem  3 

Under  the  maximum  entropy  statistical  model  induced  by  the 
decorrelation  times  as  given  in  Theorem  2  and  F' (t)  bounded  and 
continuous  in  t  rrt^,t^),  the  variational  problem  of  (20)  and 
(21)  becomes: 

-  ^1  - 

min;  J  =  dt  tr[P'V’]  (34) 

K , f , a  ^o 

where  P'  is  the  unique,  hermitian,  positive  semi-definite  solu¬ 
tion  of; 


-P'  =  l”  P'  +  P' [A'-F'-m*]  +  R'  +  I'{P'} 


(35) 


P’  (tj^)  =  0 


where 


!•  ^ 


I  0 
0  0 


(36) 


I  =  diag  , . .  . , 
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and  where,  for  any  square  matrix,  M: 


{M} 


A 


diag 


k 


(37) 


Proof 

With  the  statistical  model  of  Theorem  3,  and  the  definition 
(22)  of  stochastic  integrals,  the  results  of  Theorem  1  may  be 
used.  In  particular,  (18)  yields: 


»|)(t,T) 

=  4.'“(t+6,t) 

Ij>(t+6,T)  (j>'(t+6,t) 

a. 

\j>(t+6,T) 

=  (|.•®(T,t+6) 

R'  (t)  4)'  (T,t+6) 

b. 

where  6  >  0.  Equation  (25)  shows  that  <|)'(t+6,t)  depends  upon 
6j^(ti,t2)  only  for  tj^,t2  e(t,t+6l,  while  (|»'(T,t+6)  depends  upon 
^k^^l'^2^  only  for  tj^,t2  e[T,t+6).  Since  these  intervals  are 
disjoint  and  the  increments  of  5j^(0,t)  are  independent,  the 
ensemble  average  of  (37. a)  becomes: 

ii7(t,T)  =  E[(|)’“(t+6t)  ■ii;(t+6,T)  ())'(t+6,t)]  (39) 

where 

lJi(t,T)  =  E[»|;(t,T)] 
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Now  examine  (25).  Keeping  in  mind  that  1)  the  (J)^s  appear¬ 
ing  in  the  integrals  of  (25. c)  are  each  dependent  on  increments 
of  the  6j^(0,t)  over  mutually  disjoint  intervals,  2)  that  partial 
sums  of  (25. a)  are  almost  everywhere  convergent,  and  3)  the 
6.  (0,t)  are  Gaussian  with  variances  (33),  it  is  seen  that  the 

J>  ^ 

contribution  of  j^|2  <l'j^(t+6,t)  to  (|)'(t+6,t)  produces  terms  of 
order  6^  on  the  right  side  of  (39). 

Equation  (39)  may  thus  be  written; 


i|;(t,T)  =  E[((|)^(t+6,t) 
X  \I)(t+6,T)  ((j)^(t+6,t) 


-/^■^^dTi(j)^(t+6,Ti)F'  (T^)4,^(Tj^,t))“ 
-/^■^^dTi4.^(t+6,Ti)F'  (Tj^)  ()>;(Tj^,t))  ] 


+  0(6^) 


similarly,  using  the  expression  (25. b)  for  the  <fi^s  appearing 
above ,  we  have : 


i);(t,T)  =  E[a”  i|j(t+6,T)  A]  +  0(6^) 

A  =  I  +  A'6  +  W'(t,t+6)  +  *5W'^(t,t+6)  (Xj^) 


After  expanding  out,  rearranging  and  dividing  by  6: 


-  \  [!|)(t+6,T)  -l|»(t,T)  ] 

=  (A*  -P*  (t)  +  -^  ElW'^{t,t+6)  ])“  ■ijJ(t+6,T) 

+  ^(t+6,T){A'  -F'(t)  +  ^  ElW'^(t,t+6)  1) 

+  ^EtW'”(t,t+6)  li?(t+6,T)  W'(t,t+6)1  +  0(6) 


Next  use  Theorem  3  to  evaluate  the  above  averages,  then 
pass  to  the  limit  6+0.  Recalling  that  ^  ^(t,T)  exists  by 
Theorem  l.D,  we  obtain: 


-  ^  i(»(t,T)  =  (A'  -F'(t)  -%!•)“  'j;(t,T) 

+  l(?(t,T)(A'  -F' (t)  -hi')  +  l'{?(t,T)} 


with  I  given  as  in  (36).  Also,  (27. b)  yields  directly: 


t)  =  R'  (t) 


Finally,  integration  of  all  terms  in  (40)  over  t  e[t,tj^l  and 
use  of  (26)  gives  (35)  . 


(40.  a) 


(40. b) 


Under  the  maximum  entropy  statistical  model  we  thus  obtain 
a  modified  Lyapunov  equation  for  P'  which  must  be  appended  to 
the  variational  problem  as  a  constraint.  Clearly  as  the  decor¬ 
relation  times  approach  infinity,  the  matrix  I'  approaches  zero 
and  (35)  reduces  to  the  feuniliar  Lyapunov  equation  for  a  deter¬ 
ministic  plant.  Note  that  (35)  could  have  been  shown  by  pro¬ 
ceeding  directly  from  the  Ito  differential,  (28) .  However,  the 
method  of  the  above  proof  illustrates  the  simplicity  and  unity 
afforded  by  the  formalism  of  Theorem  1. 

2.3  Stochastic  Stability  and  the  Steady  State  Case 


Because  of  the  relative  ease  with  which  constant  gain  con¬ 
trols  may  be  implemented,  we  hence  forth  consider  only  the  steady 
state  case  and  suppose  that  K,f  and  a  are  t ime- independent .  As 
a  preliminary  step,  we  first  consider  stochastic  stability  and 
introduce  "equivalent  coefficient  matrices"  in  the  sense  of 
Kleinman  [14].  In  the  following,  the  equivalent  coefficient 
matrix  of  P'  in  (35) ,  for  o'  ample,  will  be  denoted  by: 


Ap,  [(A'  -P’  ->s-r*)”  P’+P'(A'  -F'  -hi')  + 


or  more  simply  by  A^,  whenever  clarity  permits.  By  definition, 

the  limit  of  P’  as  tj^  approaches  infinity  exists  if  and  only  if 

A—,  is  exponentially  stable.  When  this  is  so,  we  term  <,f  and 

a  admissible  controls.  Note  that  the  equivalent  coefficient 

matrix.  A?  of  the  covariance  of  X  (which  is  adjoint  to  P')  is 
U  H 

simply  A—, .  Thus  admissible  controls  imply  second  mean 

stability  (for  a  discussion  of  this  and  related  concepts  of 
stochastic  stability  see  [15]). 

More  precisely,  we  have  the  result  (see  Ref.  [16]  or  [17]). 
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Lenana  1 


SJ 


Consider: 


A®A  +  AA  +  I'{A}  +  S  ,  t 

(41) 


with  S  ^  0  (S  ^  0)  .  If  Ay^  [A^A  +  AA  +  /'{A}]  is  asymptotically 
stable. 


lim  A (t)  *  A 

^  A  “ 

tj^+OO 

where  A^  exists  as  the  unique  positive  semi -definite  (negative 
semi-definite)  solution  to: 


e 

-A  = 
A(t^) 


0  =  a"a.  +  A_A  +  J*{A.)  +  S 


(42) 


In  consequence,  under  admissible  controls,  the  steady  state 
performance: 


J,  -  lim  J/lt,-t„| 


exists  and  assumes  the  form: 
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tr  [P'V  ] 


(43) 


where  P*  is  a  positive  semi-definite  solution  of: 


0  =  lA'  -F'  -I5J']”  P'  +  P'  [A'  -F'  -55!' ] 


+  R'  +  I'{P’ } 


We  now  consider  (43)  and  (44) .  Preparatory  to  determina¬ 
tion  of  the  optimal  control  gains  we  note  the  following  results : 


Lemma  2 

Let  61...]  denote  the  first  variation  of  [...]  consequent 
upon  variations  in  tc,f  and  a  subject  to  the  admissibility  condi¬ 
tion  and  (44) .  Then  the  stationary  condition: 


6[Jg]  =  0 


determines  a  minimum  of  J  . 

s 


Proof 

The  condition  for  extremalization  of  requires  that 
6  [P' ]  vanish.  Consequently,  the  second  variation  of  P'  may 
be  computed  from  (44)  as: 


1 . 


0 


-  tA^  -P^  6^[P'l  +  6^[P')lA^  -F^  -**1’] 


+  I'{6^[PM}  + 


0 

.0 


0 

6lK]“R26[tc] 


where  A'  and  F*  denote  A'  and  F*  evaluated  with  a,K  and  £  as 
O  O  M 

determined  in  accordance  with  (45).  Since  6Ik]  R,6  1k:]  ^  0  and 

2  _  2  _  ^ 

Ag  (P']  is  asymptotically  stable,  6  [P']  is  positive  semi- 

definite  by  L^nma  1.  Since  V  s  0,  the  second  variation  of  J 

s 

is  non- negative.  □ 

In  the  following,  we  shall  concentrate  on  the  derivation 
of  the  stationary  conditions,  i.e.,  the  conditions  imposed  on 
a,K  and  £  by  the  requirements: 

6  [Jg]  =  0 

3  =  trlP'V] 

s 

0  =  (A*  -F'  -JsJ’]”  P'  +  P'  [A'  -F'  -JsJ’l 

+  R'  +  I'{P' } 


Subsequently,  at  least  for  the  case  of  full  order  compensation, 
(N^  a  2n) ,  the  conditions  under  which  A^,  is  asymptotically 
stable  will  be  established. 
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3.  THE  MEAN-SQUARE  OPTIMAL  FULL-ORDER  COMPENSATOR 

3.1  Derivation  of  the  Stationary  Conditions 

Here  and  in  the  remainder  of  this  report  we  consider  the 
case  of  full-order  dynamic  compensation  for  which  N^  =  2n  in 
(46) .  Leaving  aside  the  question  of  admissibility  for  the 
moment,  let  us  now  derive  the  stationary  condition  (4 6. a). 

To  handle  the  constraint  imposed  by  (46. c)  most  expedi¬ 
tiously,  we  introduce  the  hermitian  multiplier  matrix, 

Q'  e  and  form  the  Hamiltonian: 


H  =  tr[P'V'  +  ([A'  -P'  ->5l']^  P'  +  P' [A'  -F'  -JsJ'l 

+  R'  +  I'{P'  })Q'  ] 


(47) 


Then  (46. a)  reduces  to  the  requirement  that  the  first  variation 
of  H  consequent  upon  unrestricted  variation  in  a,K,f  and  P' 
vanish.  Partitioning  P'  thus: 


P* 


(48) 


with  similar  notation  for  Q' ,  we  obtain  the  following  specific 
conditions : 
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on 

ap. 

S 

[A'-F'-y  IQ-  +  Q'  [A' 

0 

-F'-Jjl'  ]  +  J*  {Q'  }  +  V 

a. 

dH 

dK 

+  ° 

b. 

dH 

5f 

s 

+  PqfVj]  .  0 

c. 

dH 

da 

= 

0 

d. 

From  (49. a),  Q'  can  now  be  identified  as  the  covariance  of 
the  augmented  state.  Furthermore,  it  is  easily  checked  that  by 
virtue  of  (46. c)  and  (49.a,b,c),  (49. d)  is  an  identity.  This 
manifests  the  well-known  result  that  the  stationary  conditions 
furnish  no  determination  of  the  compensator  dynamic  matrix,  a. 
Indeed  it  is  fortunate  that  we  may  choose  a  so  as  to  simplify 
(46. c)  and  (49)  very  greatly. 

In  preparation  for  the  main  result,  define  the  state  trans' 
formation: 


( 


so  that  q  is  the  observation  error  if  we  view  the  compensator 
as  a  full-order  observer.  Noting  that 
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pi  _  T  P  T 


Q  t  =  T  Q 


partition  P,  and  Q  in  accordance  with  (48) ; 


P 


Finally,  with  the  notation: 


y  -hi 


we  have: 

Theorem  4 


With  the  above  definitions  and  the  choice: 


the  specifications: 


31 


B 


and 


K  =  R-^  6»  Pj 


f  =  5,rV 


b. 


■( 


(55) 


0  = 


(W  -Bic)”  Pp  +  Pr<Vi«-Bic)  +  r,  +  k®R,<  +  i{Pp+  ?  } 
m  1  i  ^  q 


0  = 


(56) 


+  J{Q^}  + 


0 

0 


(y^-fY)“  Pq  +  Pq(V^^>  •*■ 


(57) 


(y„-fY)Qq  +  Qq(y^-fY)”  +  liQ^}  +  +  fv^f 


H 


identically  satisfy  the  stationary  conditions  (46. c)  and  (49). 
The  proof  is  given  in  Appendix  1. 

with  the  choice  (53) ,  system  (17)  assumes  the  form: 


C  =  y?  -Bicq  + 


4  *  y^q  “BKq  +  fY(?-q)  +  fWj 


b. 


■( 


(58) 


32 


■lAa . .  Xii  .Vi  v'n  >. 


so  that  the  compensator  structure  is  as  shown  in  Figure  1. 
Furthermore,  referring  to  (28) ,  the  mean  value  of  C  is  seen 
to  satisfy: 

t  =  I  -3<q 

Thus,  it  is  clear  from  (58. b)  that  with  a  chosen  as  in  (53) ,  thi 
compensator  is  a  full-order  observer  of  the  mean  state. 

It  is  also  evident  that  under  (53),  Eqs.  (55),  (56)  and  (5 
comprise  the  only  essential  stationary  conditions.  Introducing 


Corollary  1 

With  the  full-order  compensator  form  given  in  (53.b),  the 
gains : 


K  =  R2^  b”  P 


x:  «  H  -1 

f  =  QY  V2 


/V  /V 


where  P,  Q,  P,  Q  2  0  satisfy; 


-H 


0  =  y"  P  +  P  u  +  J{P  +  P}  +  r.  -P  a  P 

mm  1 


ttH 


“  =  V„  Q  +  Q  +  X{Q  +  Q)  +  -Q  o  Q 


0  =  (y^jj  -Q  a)“  P  +  P  (y^j^  -Q  a)  +  P  a  P 


/\  At 


H 


0  =  -a  p)  Q  +  Q  (y^^  -a  p)  +  Q  a  Q 


(6 


a. 


(6 


b. 


a. 


(e 


b. 


yield  an  extremum  for  the  steady-state  performance,  J  . 

9 

The  particularly  simple  form  of  (62)  and  (63)  is  to  be 
noted.  Observe  that  in  the  deterministic  plant  case  (J  =  0) , 
(62)  form  the  only  restrictions  and  reduce  to  the  fcuniliar 
uncoupled  regulator  and  observer  Riccati  equations.  Also  in 

^  A 

the  absence  of  the  terms  J{P}  and  J{Q},  (62)  are  of  the  form 
of  the  stochastic  Riccati  equation  arising  from  the  regulator 
problem  treated  earlier.  In  general,  however,  (62. a)  and  (62.b) 


are  coupled  through  J{P}  and  j{Q}  so  that  the  effects  of  modal 

frequency  uncertainties  preclude  the  separation  principle.  Also 

the  coupling  terms  effectively  augment  the  diagonal  elements  of 

r,  and  v. ,  thereby  demanding  greater  control  authority  and 

stabilization.  Moreover,  from  (63),  p  and  Q  tend  to  increase 

H  H 

with  the  controller  input  <  R2K  and  its  dual  f  V2f. 

We  may  say  that  P  and  Q  represent  error  "leaking  through" 
the  regulator  to  the  observer  and  vice-versa  by  virtue  of  pareun- 
eter  uncertainty. 

To  summarize,  we  consolidate  the  above  results  and  simplify 
the  requisite  conditions  for  admissibility  of  the  controls; 

Theorem  5 

Suppose  that  P,  Q,  P,  Q  exist  as  positive  semi-definite 
solutions  to  (62)  and  (63).  Then  with  constant  gains  given  by 
(61)  ; 

A.  (il  -fy)  ,  (y  -fy)  »  (ir  “6k)  and 

Rt  tCl 

(y  -3k)  are  asymptotically  stable. 

B.  The  stochastic  system  (58)  is  second  mean  exponen¬ 
tially  stable  and  almost  surely  exponentially  stable, 
i.e.,  the  control  is  admissible. 

C.  Eq.  (61)  yield  a  minimum  of  the  steady  state  per¬ 
formance,  7  . 

8 

The  proof  is  given  in  Appendix  2. 


At  this  point,  it  is  well  to  recapitulate  the  above  develop¬ 
ments,  indicating  where  the  derivation  is  subject  to  significant 
generalization.  First,  we  chose  to  treat  the  problem  of  mean- 
square  optimal  dynamic  compensation  of  structural  systems  with 
uncertainties  in  the  open  loop  frequencies  within  the  dual  form¬ 
ulation,  involving  the  augmented  co-state  or  expected  cost  matrix. 
This  gives  rise  to  the  statement  embodied  in  Eqs.  (20)  and  (21) 
which  is  applicable  to  any  statistical  model  of  parcimeter  uncer¬ 
tainties.  Next  we  imposed  the  minimum  information  statistical 
model  of  frequency  uncertainties  to  obtain  the  variational  prob¬ 
lem  of  Eq.  (46)  which  involves  a  variational  constraint  imposed 
by  the  stochastic  Lyapunov  equation  for  the  expected  cost.  Treat¬ 
ment  of  less  restricted  types  of  pareuneter  uncertainty  would  pro¬ 
ceed  analogously  to  obtain  an  appropriately  generalized  stochas¬ 
tic  Lyapunov  equation. 

Assuming  admissibility  of  the  controls,  and  specializing  to 
the  steady  state,  full-order  case,  we  derived  the  optimality  con¬ 
ditions  associated  with  the  variational  problem  of  Eqs.  (46) . 
Choice  of  the  compensator  as  a  full-order  observer  for  the  mean 
state  is  found  to  permit  great  simplification  of  these  conditions 
resulting  in  Eqs.  (61)  through  (63).  Clearly,  the  case  of  time- 

varying  controls  may  be  handled  in  a  similar  manner.  Moreover, 

N 

we  may  anticipate  that  assuming  q  eCQ,  <2n  and  proceeding 
as  in  this  section,  the  optimality  conditions  for  reduced-order 
dyneunic  compensation  may  be  derived.  Once  again,  judicious 
choice  of  the  compensator  dynamic  matrix,  a,  will  result  in 
drastic  simplification  and  yield  appropriately  modified  versions 
of  Eqs.  (61),  (62)  and  (63). 


Finally,  Theorem  5  shows  that  existence  of  positive  semi- 
definite  solutions  to  the  optimality  conditions,  Eqs.  (62)  and 
(63) ,  guarantees  the  minimum  property  and  second  mean  and  almost 


siire  exponential  stability  for  the  closed  loop  system.  Thus  the 
new  approach  to  full-order  compensation  holds  the  promise  of 
ensuring  robust-  stability  under  a  design-conservative  statisti¬ 
cal  model  of  parameter  uncertainties.  The  issues  of  existence 
and  uniqueness  are  immediately  addressed  in  the  following  section. 

3.2  Existence  and  Uniqueness  of  Solutions  to  the 
Optimality  Conditions 

The  results  of  this  section  represent  an  extension  of  the 
earlier  work  on  full  state  feedback  regulation  to  the  more  gen¬ 
eral  problon  posed  by  Eqs.  (62)  and  (63) .  As  in  the  classic 
work  of  Wonham  [17,18],  or  the  developments  of  Merrieun  [16],  the 
basic  technique  is  to  establish  the  existence  of  bounded  mono¬ 
tone  sequences  of  positive  semi -definite  hermitian  maps.  In 
this  connection  the  following  result  is  essential  [19] : 

Lemma  3 

Every  sequence,  {Xj^},  of  hermitian  positive  semi-definite 
matrices  bounded  below  (above,  resp.)  with  "^i^  negative 

semi -definite  (positive  semi-definite,  resp.)  for  each  i  con¬ 
verges  to  a  positive  semi-definite  limit. 

In  addition,  we  shall  need  a  preliminary  lemma  on  the 
familiar  Lyapunov  equation  [16,20]. 

Lemma  4 

Consider : 


A**  P  +  P  A  +  S 


0 


(64) 


A.  If  S  ^  0  0,  resp.)  and  A  is  stable,  (64)  has  a 

unique  solution  P,  and  P  ^  0  0,  resp.).  If,  in 

if 

addition  (S  ,  A)  is  reconstructible,  P  >  0  ( <  0,  resp.) 

B.  Suppose  P  2  0  0,  resp.) ,  S  >  0  0,  resp.)  satisfy 

(64)  and  (S**,  A)  is  detectable.  Then  A  is  stable. 

With  regard  to  stochastic  Lyapunov  equations  of  the  form: 


0  =  a“p  +  P  a  +  J{P}  +  S 


(6: 


we  can  state  [14] : 

Lemma  5 

If  Ap  [A*^  P  +  PA+J{P}]  is  asymptotically  stable  and  S  ^  0 
(50,  resp.)  then  the  unique  solution  of  (65)  is  positive  (nega¬ 
tive,  resp.)  semi-definite. 

Finally,  we  must  recapitulate  earlier  results  for  the 
stochastic  Riccati  equation  of  the  regulator  problem  [7] : 

Lemma  6 

Consider: 


tl 

=  A  P  +  P 
m 


+  I{p}  +  r  -P 


a  P 


0 


where  r  ^  0  and  I  and  a  are  as  defined  by  (36.b)  and  (60. a). 
Under  the  conditions: 


(A,  a  )  stabilizable 

(r^.  A)  detectable 


(66)  possesses  a  unique  positive  semi-definite  solution  P,  and 
(A  -aP)  /  (A  -hi  -aP)  and  AplA^  P  +  P  A^  +  r{P}]  are  asymptoti¬ 
cally  stable.  Moreover,  suppose  that  P^^  and  the  solu¬ 

tions  of  (66)  with  r  =  r^^  and  r  =  r2,  respectively.  Then 

^1  ^  ^2  ^  ^2* 

The  dualization  of  Lemmas  4,  5  and  6  is  immediate  and  need 
not  be  considered  explicitly. 

To  simplify  this  initial  development,  attention  is  hence¬ 
forth  restricted  to  the  case: 

—  h  —  h 

(U/  o^)  ,  (vi,  )  controllable 

h  —  —h  — 

(r^^^,  v)  ,  (a^,  y)  reconstructible 


although,  as  subsequent  nixmerical  results  suggest,  it  is  likely 
that  these  conditions  can  be  weakened  considered)ly. 

Within  the  above  restrictions,  consider  the  sequences 

A 

Pi#  Pj^  and  defined  by; 


40 


0 

0 

0 

0 


u**P.  +  P.u  +  r,  + 
mi  1  m  1 


^mQi  +  + 


J{P^} 

liQ^} 


-Pi  a  Pi  +  i{Pi} 

-Qi  a  Qi  +  J{Qi} 


^i+l^^m  -Q^a)”  P.^i  +  P,  a  P 


m  1 


i+1 


Q.^,(y  -aP,)®  +  (y  -aP.)  Q_,  +  Q.  a  Q. 

1+1  m  1  m  1  '^1+1  '^1  '^1 


a. 

b. 


a. 


b. 


I  (67: 


I  (68: 


for  all  i  ^  0.  The  following  results  on  the  boundedness  of  these 
sequences  is  essential. 

Theorem  6 

Under  the  conditions  stated  above,  there  exist  bounded, 
positive  definite  P^,  P^,  p'^,  P^,  Q^,  Q*”,  q'^,  with 


,u 


i  P 


u 


''  a 


^u 


q“  ^  q'^ 


such  that: 


(69: 


in  conjunction  with  (67)  and  (68)  implies 


■  ♦..fWlHf 


(y  -aP.)  ,  (y  -Q.  a)  stable 
mi  mi 


and  for  i  and 


<  P, 


s  Q 


u 


for  1^1. 

The  proof  is  contained  in  Appendix  3. 

This  prepares  the  way  for  the  main  conclusion: 

Theorem  7 

With  y^,  I,  a  and  o  as  defined  previously  and  the 

conditions : 


(M/ 

,  (y,  Vj^S 

controllable 

y) 

,  {a^,  y) 

reconstructible 

Eqs.  (62)  and  (63)  possess  unique  positive  definite  solutions 

for  P,  P,  Q  and  Q.  Moreover,  (iL  -crP)  and  (11  -Qa)  are 

m  m 

asymptotically  stable. 


Proof 


First,  let  us  stipulate  that  P  and  Q  be  chosen  in  accord 

O  O  /s 

emce  with  (69)  and  Theorem  6.  Then,  for  all  i  s  0,  P . ,  P. ,  Q. 

/\  i  i  i 

and  as  defined  by  (67)  and  (68)  are  positive  definite  and 

(y  -aP.)  and  (y  -Q.o)  are  asymptotically  stable, 
m  1  m  1 

Defining: 


Pi 

4 

As,  /\ 

-P. 

1+1  1 

Qi 

A 

As  A 

Qi+1  -Qi 

and 


Pi 


P.^,  -P. 
1+1  1 


0^1  -°i 


manipulation  of  (67)  and  (68)  yields 


0 


■  <%  -Oi  -Qi  ® 

1  1  1-1  i-1 


(72. a) 


■  ^p  +  ^p  O'  P.  1 

1-1  ''i-1  ^  i-1  *^i-l  ^“•‘' 


*i-l 


0  = 


\''‘«  *  '%  -<*i)  "q.  -!:q.  /  , 

^  1  ''i-l  “i-l 

>'•  A 

-  [Q^Sp  a  +  a  I  Qil  +  ^  Q.-  +  Q.a  Z 

^  1-1  ^i-1  ^  ^i-l  ^  ^  Qi-1 


(72. b) 


and 


'“m  ^P^  +  '%  -'^i) 


+  l{Ep  +  Z_  }  -E_  a  E- 
^i  ^i  ^i  ^i 


a. 


-.H 


‘“m  -0l  \  *  %  (*■„  -Qi  o) 


"  ^«Qi  *  ^Qi>  -^Q/ 


b. 


(73) 


Now  suppose  Zp  i  0  and  Z^  so.  Considering  (7 3. a), 
it  is  clear  that  siici  y  -oP.  is^siable,  {Z  -aP.  ,,a^)  and 

J)  *  ”^^i-l^  and  detectable, 

respectively.  Then,  by  Lemma  6,  Ep  exists  as  the  unique 
positive  semi-definite  solution  of  (73. a).  Similarly,  from 
(73.b),  Lemma  6,  and  the  assumption  Z^  s  0;  E  is  negative 
semi-definite. 


44 


Since  -Q.  ^ 


^  0  and  a  ^  0,  use  of  (68. a)  and  Lenuna  4 


shows  that  [P.  Z_  a  +  a  P.]  is  negative  semi-definite. 

1  Qi_l  1 

Similarly  [P .  ,  a  E  a  P.  ^  0  so  that  the  second 

^i-1  ^i-1 

line  of  (72. a)  is  positive  semi-definite.  Consequently,  since 
(il  -Q^a)  is  stable,  Zp  is  the  unique  positive  semi-definite 
solution  of  (72. a)  by  virtue  of  Lemma  4. A.  Analogous  reasoning 


on  (72. b)  shows  that  Z_  s  0. 

^i 

Thus,  on  the  assumption  that  Zp  ^0  and  Zq  ^  0,  we 

^i-1  ^i-1 


have  shown  that  Z  ^  0,  ^  0,  Z 

^i-1  ^i-1  ^ 

Induction  on  i  shows ; 


^  0  and  Z 


i-1 


Qi-1 


s  0. 


1 

1 

IV 

o 

o 

VI 

s  0 


and  Q.  -Q  ^  0. 
1  o 


for  all  i  provided  that  P^^  -P^ 
choice  (69)  this  must  be  so,  for  by  Theorem  6,  P 

/V  /s  y\ 

Ql  ^  Q  =  Qq. 


Under  the 

/N  j  ^ 

P  =  P  and 
o 


Thus  (67)  ,  (68)  and  (69)  define  Pj^  and  P^  as  positive  defi¬ 
nite  nondecreasing  sequences  and  Q.  and  as  positive  definite 
non- increasing  sequences.  Moreover,  by  Theorem  6,  these  sequences 
are  bounded  both  from  above  and  from  below.  Therefore  Lemma  3 
implies  that  the  sequences  P^,  P^,  and  possess  positive 
definite  limits,  which  by  virtue  of  (67)  and  (68)  satisfy  (62) 
and  (63) . 


45 


Furthermore,  [u  -oP.)  and  (y  -Q.cr)  are  asymptotically 
mi  mi 

stable  so  that  lim  (y  -CP^)  and  lim  (y  -Q.o)  are  also  stable. 

i+co  .  ^  .  i+co  ^ 

Finally,  since  lim  P.  and  lim  Q.  exist  the  results  of  Lemma  6 

i+oo  ^  i+oo  ^  ^ 

immediately  apply  to  (67) .  Thus  buth  y  -aP  and  y  -Qa  are 
asymptotically  stable. 

To  establish  uniqueness,  suppose  that  (Pj^,  Pj^,  and 

(P2,  '^2'  ^2'  sets  of  positive  definite  solutions  of 

(62)  and  (63).  With  =  Pj^  -P2  and  Zp  =  -P2,  (62.  a)  may  be 

manipulated  to  yield: 


0  -  (y„  -crPj)"  Ep  +  Ep(y„  -OP^)  .  /{Ep) 


+  j(Zp}  -Ep0  Ep 


0  -  (y„  -OPj)”  Ep  +  Ep(M^  -0Pj)  .  KEp) 


+  r{Zp)  +  EpO  Ep 


Employing  the  same  reasoning  as  used  in  the  Appendix  for 
deriving  bounds  on  P,  it  is  readily  seen  that  J{Zp}  <  IpC  £p. 
Thus,  since  in  (74. a)  is  asymptotically  stable,  Ep  is  nega¬ 
tive  seml-def iSite  by  (74. a)  and  Lemma  5.  Likewise,  (74. b) 
implies  Ep  2  0.  Since  Ep  is  both  positive  and  negative  semi- 
definite,  we  conclude  that  Ep  =  0.  Similarly,  =  ©2.  Finally 
(P.  -Po)  and  (Q.  -Qo)  found  to  satisfy  the  homogeneous  forms 

of  (63),  whence  P,  =  P-  and  Q,  =  Q-.  □ 


As  is  clear  from  the  above  proof,  (67) ,  (68)  and  (69)  give 
a  computational  procedure  for  the  numerical  solution  of  (62)  and 

'•0  ''ll 

(63) .  When  P  and  Q  are  evaluated  in  accordance  with  Eqs. 

(A. 3.1)  through  (A. 3. 4)  of  the  Appendix,  monotone  convergence  is 
guaranteed  (with  P^  non-decreasing  and  Q^,  non-increasing) 
It  is  interesting  to  note  that  if  the  conditions: 


=  P‘ 


Qo  =  Q 


are  taken  in  place  of  (69) ,  slight  modification  of  the  above 
proof  again  shows  convergence  of  (67)  and  (68),  but  with  P . ,  P. 

A 

non-increasing  and  Q^,  non-decreasing. 

Theorem  5  is  now  directly  applicable.  Thus  it  is  possible 
to  summarize  definite  conclusions  with  regard  to  stochastic 
stability  as  follows: 

Corollary  2 

A  /V 

Under  the  conditions  of  Theorem  7,  let  P,  Q,  P  and  Q  be  the 
unique  positive  definite  solutions  of  (62)  and  (63) .  Then  with 
the  constant  gains  given  by  (61) : 

A.  (TIjj  -3k)/  (II  -3k),  (Ujjj  -fy)  and  (II  -fy)  are  asymp¬ 
totically  stable. 

B.  The  stochastic  system  (58)  is  second  mean  exponen¬ 
tially  stable  and  almost  surely  exponentially  stable. 

C.  ^  (61)  minimizes  the  steady  state  performance  index 

defined  by  (43) . 


Thus,  under  the  minimum  information  model  of  modal  frequency 
uncertainties,  the  compensator  design  approach  embodied  in  Eqs. 
(61)  through  (63)  minimizes  a  measure  of  performance  defined  over 
the  entire  pareuneter  ensemble  and  guarantees  closed  loop  stochas¬ 
tic  stability. 

3.3  Asymptotic  Properties  for  Large  Uncertainties 

The  previous  section  established  the  existence  and  unique¬ 
ness  of  solutions  to  (62)  and  (63)  under  the  controlability  and 
reconstructability  conditions  of  Theorem  7  and  for  all  positive 
I.  It  is  thus  natural  to  inquire  what  behavior  P  and  Q  attain 
for  large  uncertainties,  i.e.,  for  very  small  decorrelation  times 
This  case  represents  the  situation  in  which  very  great  a  priori 
uncertainty  exists  regarding  the  values  of  all  modal  frequencies. 

In  connection  with  the  stochastic  Riccati  equation.  Theorem 
16  of  Ref.  [17]  immediately  yields  the  following: 


Lemma  7 

Assume  the  conditions  of  Theoran  7.  Let  P  and  Q  be  the 
unique  positive  definite  solutions  of  (62)  given  any  positive 
semi-definite  P  and  Q.  Introduce  a  positive  scaling  parameter, 
J.  into  J: 


I 


lim  P 
J  +  oo 

lim  Q 
Jf® 


where  the  diagonal  matrices  P*  and  Q*  are  uniquely  determined  by: 


0  =  p”  P*  +  P*  y  +  I{P}  +  {r^}  -  {P*}^  {a} 
0  =  y  Q*  +  Q*  y“  +  I{Q}  +  {v^^}  -  {Q*}^  {a} 


This  result  readily  leads  to  the  main  conclusion 
Theorem  8 

Let  P  and  Q  be  the  unique  positive  definite  solutions  of 
(62)  and  (63)  under  the  conditions  of  Theorem  7  and  define  I  as 
in  (75).  Then: 


lim  P  =  P*  ^  (-2  Rey)"^  {r,} 

J+oo 

lim  Q  =  Q*  =  (-2  Rey)“^  {v, } 


Proof 

By  Theorem  7,  P  and  Q  exist  as  positive  definite  solutions 
of  (62)  and  (63)  for  all  1^0  and  therefore  possess  bounded 
limits  P*  and  Q* ,  respectively,  as  J  increases  without  bound. 
Then  (77)  yields: 


49 


0 

0 

—  h 

Furthermore,  since  J  >  0,  ((Pa  P)  ,J)  is  reconstructible 
— 

and  (J, (Q  a  Q)  )  is  controllable  for  all  P,  Q  >  0.  Then  (63. a) 
gives : 


=  y®  P*  +  P*  y  +  t{P*}  +  {rj^}  -P*^  {a} 
=  y  Q*  +  Q*  y**  +  J{Q*}  +  {v,  }  -Q*^  {a} 


(79) 


{P*}  4  lim  {P} 

J  +  oo 


l“^  lim  {P  o  P} 

Jfoo 

{p*}^  {a} 


(80. a) 


Similarly,  from  (63.b) : 


{Q*}  ^  lim  {Q} 

J+oo 


{Q*}^  {a} 


(80. b) 


and  (78)  follows  by  substitution  of  (80)  into  (79) .  □ 

Thus,  as  modal  frequency  uncertainties  increase  without 
bound,  the  optimality  conditions  possess  very  simple  asymptotic 
solutions  given  by  closed  analytical  expressions.  Referring  to 
Theoron  16  of  Ref.  [7]  it  must  be  noted  that  for  the  full  state 
feedback  regulation  problem,  diagonalization  of  the  expected 
cost  matrix  gives  rise  to  a  rate  feedback  control  law  which  is 
inherently  stable  for  all  values  of  modal  frequencies.  Thus 
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implies  that  only  the  estimated  modal  valocities  are  used  as 
control  input. 


In  this  case,  it  is  also  of  interest  to  investigate  more 
particularly  the  asymptotic  form  of  the  control  law  for  large  J. 
For  this  purpose,  suppose  that  y,  g,  and  y  are  the  modal  par- 
cuneters  for  a  particular  realization  of  the  structure  as  distinct 
from  the  nominal  or  design  values,  y,  3  and  y.  Then  the  system 
equations  for  the  actual  structure  in  connection  with  the 
designed  compensator  are: 


i  =  y  5  -6  K  q 
q  =  (v^  -gK  -fy)  q  +  f  Y  5 

Now  as  J  >  0  increases,  the  compensator  poles  recede  into  the 
left  half  plane  owing  to  the  term  -hi  in  iljj*  Thus,  for  large  J 
we  may  consider  a  singular  perturbation  expansion  for  q  as  a 
functional  of  In  the  first  approximation  -  i.e.,  to  first 
order  in 


q  ~  21"^  f*  Y  5  +  0(1“^) 


(82 


where  f*  =  lim  f.  Then,  to  first  order,  in  singular  perturba- 

J  A  QQ 

tions,  the  controlled  system  is  given  by: 


y  5  -2  a(P*  I  *  Q*)  a  C 


~  A  ~  -1  H 

a  =  6  Rj  6 


-  A  H  „-l  ~ 
a  ^  y  ^2  ^ 


from  (82)  and  (81. a). 

To  zeroth  order,  (83)  shows  that  large  uncertainties  tend 
to  remove  the  control  altogether.  However,  to  first  order  we 
have  the  following  property: 

Lemma  8 

Consider  (83)  with  P*  and  Q*  as  given  by  (78).  If  y  = 

Y  *  8  ,  and  v.  =  vH.  (where  v  is  a  positive  scalar)  then  (83)  is 
asymptotically  stable  for  all  y  with  Key  <  0  and  all  Br  and  the 
control  is  a  rate  output  feedback  law. 


Proof 


Let  £  =  C**  C-  Then,  from  (83)  : 


2^"  [ReC  -n  K 


T  »  a(P*  l"^Q*)  0  +  c  (P*  J“^Q*) 


In  the  case  y  =  y  /  Y  =  p  and  V2  ■  VR2  we  have  a  *  va  so 
that  r  is  hermitian.  Further,  since  P*  J~^Q*  is  positive  and 
diagonal,  F  2  0.  Then,  with  Rey  <0,  £  <0  for  all  ^  ^  0.  Thus 
£  is  a  Lyapunov  function.  That  this  case  represents  a  rate 


output  feedback  law  is  easily  verified  by  transforming  (83)  back 
to  the  modal  coordinate  basis  in  accordance  with  relations  (7) 
and  (87) .  □ 

In  other  words,  for  large  modal  frequency  uncertainties  and 

^  H 

Co-located  rate  sensors  (y  =  3  »  Y  =  8  ) »  the  asymptotic  control 
is  a  rate  feedback  law  proposed  by  Balas  [4]  which  is  stable  in 
the  presence  of  errors  in  all  modal  parcuneters. 

Analogous  results  are  to  be  expected  when  uncertainties  in 
low  frequency  modes  are  small  while  modelling  accuracy  degener¬ 
ates  for  modes  of  increasing  order.  In  this  instance  we  antici¬ 
pate  that  the  control  designed  according  to  Eqs.  (61)  through 
(63)  will  necessarily  approach  the  asymptotic  form  for  the  high 
frequency,  poorly  known  modes,  while  resembling  the  determinis¬ 
tic  plant  compensator  design  for  the  low  order  modes. 
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4. 


FULL-ORDER  COMPENSATION:  NUMERICAL  PROCEDURES  &  RESULTS 


4. 1  Computational  Methods 

For  computational  purposes  it  is  advantageous  to  work  with 
(61)  through  (63)  in  the  original  modal  coordinate  basis  rather 
than  in  the  complex  form  given  above.  Letting: 


p 

A 

*-lH„.-l 

4  P4 

,  P 

A 

4  P4 

s 

4 

4  0  4® 

,  Q 

4 

4  Q  4" 

(85) 


2Uld 


~  A  ;.-l 
q  =  4  q 


K  =  k4 

F  =  4  f 


-1 


(86) 


syst^  (17)  assumes  the  form 


X  =  (A  -«■  a  ( t) )  X  -BKq  +  Wj^ 
q  =  (A  -JjDq  -BKq  +  F  C(x  -q)  +  F  (#2 


(87) 


Using  (5)  through  (8),  (14. c)  and  (85),  relations  (61) 
through  (63)  become  (we  shall  suppress  the  tildes  of  P,  Q,  etc. 
in  (85))  : 


F 


Q  C*  V2 


0  =  P  P 

»  - 

0  =  (A.  -Q  E)'^  P 

m 

0  =  (S  -E  P)  Q  + 

m 

where : 


A„  ^  A  -HI 
m 


—  1  T 
B  R2^ 


c*^  v"^  c 


Dp  [I,  p  +  p]  +  -: 

Dq[I,  Q  +  Q]  +  -( 


+  P(A  -Q  E)  +  P  E  P 
m 

Q(A,  -E  P)  +  Q  E  Q 
m 


and  where  for  any  square  matrix,  M: 


(9 


Dp[J,Ml 


Dq[I,M1 


s  block-dlag 
Ic"  X  /  •  •  •  ^ 


»  block-diag 

k"lf • • • /R 


V2k,2k^ 

Lo 

.)j 

^'^«2k-l,2k-l 

®*2k,2k^ 

1 

0 

I ) 

We  also  note/  for  future  reference,  that  the  asymptotic 
solution,  (78) ,  for  large  uncertainty  levels  becomes: 


P* 


Q* 


block-dlag 
k»l, . . .  ,n 


block-diag 
k™l , • • • ,  n 


^Vk  2k-l,2k-l  •^2k,2k 

— ] 
^Vk  2k-l,2k-l  ■^2k,2k 


1 

0 


(9; 


in  the  modal  coordinate  basis. 

With  the  above  expressions,  the  discussion  of  computational 
techniques  may  proceed.  First,  note  that  the  proof  of  Theorem  7 
gives  rise  to  a  convergent  sequence  of  approximations  to  P,  Q, 

A  /s. 

P,  Q  and  yields  the  following  algorithm: 

Theorem  9 

Under  the  conditions  of  Theorem  7  and  denoting  the  positive 

As  As 

definite  solutions  of  (89)  and  (90)  by  P,  P,  Q  and  Q: 


(P,  P,  Q,  Q) 


lim  (P.  ,  P.  ,  Qj^,  Qj^)  = 
k+“ 


where  the  sequences  {P]^}»  ^®k^  ^®k^  defined 

( for  k  =  0 ,...,“)  by ; 


0 

0 


«k  ^  5^  +  Vj  +  DgU, 


"■k  +  ''k> 
°k  *  °k' 


0  = 


vxT 


-  '’k+i  -Ok  ^  <^„  -Ok  '■k+i 


°  -  Ok^i  (%  -2  Pfc)  +  (\  P^)  Q^^.1  +  Qk  0^ 


with  either; 


or 


Po  =  '■  '  Oo  -  Q“ 


Pq  =  p""  p  Qq  = 


''a  ''t  ^u  ''u 

Furthermore,  P  ,  Q  ,  P  and  Q  are  defined  as  the  positive 

definite  solutions  of 


0 


a. 


0  -  sj  pU  ^  pU  ^  ^  pU  ^  ;uj  _pU  j;  pU 

0  »  ■*■  %  °'^  ■*■  ^^  ^  °* 
0  =  q"  +  Q^  ^  +  Vi  +  Dq[j,  q"  +  5^1  -q'^  z  q’^  d. 

0  =  “Q^  “Q^  X)"^  P^  +  P^  Z  P^  a. 

0  =  P^  +  p""  A^  +  R^  +  Dp[J,  P^]  -P^  E  P^  b. 

0  =  5Na^-zp")^+  + 

0  =  ^  q'  +  q'  ^  +  ^1  +  °pt^'  q'i  -q"  ^  q'  <'• 


The  proof  is  inunediate  from  that  of  Theorem  7  and  use  of 
the  transformations  (7),  (8),  (85)  and  (86). 

The  above  computational  scheme  requires  at  each  step,  the 
solution  of  Lyapunov  equations,  standard  Rlccati  equations  and 
stochastic  Riccati  equations  of  the  form  considered  in  [7]. 
This  is  obviously  the  case  for  (95) ,  (96)  and  (100) .  To  see 
that,  at  most,  solution  of  stochastic  Riccati  equations  is  re¬ 
quired  for  computation  of  the  starting  values,  consider  (99) . 

Setting  P  ^  P**  +  p'^,  addition  of  (99. a)  and  (99.b)  produces 
s 

the  stochastic  Riccati  equation: 


(99) 


(100) 


0 


(101) 


P  A  +  P  +  R,  +  D„[J,  P  ] 
s  m  ms  1  P's 


This  defines  P  >  0  uniquely,  and  (99.b)  may  be  written 
s 


Thus,  p'^  is  the  unique  positive  definite  solution  of  an  ordinary 
Riccati  equation.  Finally,  solution  of  the  Lyapunov  equation 

/V 

(99. a)  yields  P^.  An  analogous  scheme  may  be  used  for  determi¬ 
nation  of  and  from  (99.c,d). 

Although  convergence  of  the  sequences  defined  by  (95) 
through  (100)  is  assured,  the  numerical  procedure  involves 
stochastic  Riccati  equations  at  each  step,  and  these,  in  turn, 
dCTiand  an  iterative  method  of  solution. 

In  place  of  the  above  rather  cumbersome  method  of  solution 
we  may  introduce  a  much  more  direct  Iterative  method  based  upon 
the  following  computational  sequence  (k  =  0,  1,  2,...): 


0 


0 


'■k>  Vl  ^  Vl<\,  "k'  ^  ’>1 


+  DpIJ,  +  P^J  +  E  P„ 


(103) 


V«kil  *  Ok+l<^n  V''  *  °k  ^  °k 


" 


-Ok  '^’Qk+l  *  O^^iCJV  -Q^  £)  *  v^ 


+  DqU,  Q,^  +  Q^l  *  L  Qy. 


°  ’  <%  -Ok  Vl  *  -Ok  *  "k  '’k  »• 


(104) 


with  starting  values  defined  by; 


P„  =  Q„  =  0 

o  o 


A*  +  P  A  +  R,  -P  S  P^  =  0 

o  o  loo 


A  Q„  +  Q^A  +  V,  -Q^  z  Q„  =  0 
O  O  1  o  o 


(105) 


(106) 


Note  that  (103)  and  (104)  entail  only  the  solution  of 
Lyapunov  equations.  Because  of  its  evident  convenience,  this 
approach,  in  preference  to  that  of  Theoren  9,  has  been  imple> 
mented  computationally.  Although  proof  of  convergence  r^ains 
the  object  of  investigation,  the  sequence  defined  by  (103) 
through  (106)  have  been  found  to  be  convergent  in  all  numerical 
studies  performed  to  date. 

In  the  specific  implementation  of  (103)  through  (106)  used 
to  obtain  the  numerical  results  discussed  in  the  following 
sections,  the  iterative  sequence  is  terminated  once  all  diago¬ 
nal  entries  of  P  and  Q  converge  to  within  a  given  tolerance  - 
i.e.,  given  e  >  0,  solution  of  (103)  and  (104)  is  carried  up  to 
k  ■  ( e) ,  where 


<■  • '  •  •  •  ’  • 
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The  Lyapunov  equations  (103)  and  (104)  are  solved  by  the 
real  Schur  vector  approach  following  Kitagawa  [21]  while  the 
Riccati  equations  (106)  are  solved  by  the  real  Schur  vector 
technique  of  Laub  [22] . 

The  major  computational  burden  in  solving  Lyapunov-type 
equations  is  in  reducing  the  stability  matrix  to  real  Schur 
forms.  Hence,  computation  is  reduced  by  half  when  only  one 
reduction  is  used  to  solve  (103. a)  and  (I03.b)  and  one  for 
(104. a)  and  (104. b).  Clearly,  the  solution  sequence  of  (103) 
and  (104)  also  reduces  storage  requirements. 

4. 2  General  Description  of  Design  Studies 

In  the  remainder  of  this  report,  we  discuss  application  of 
the  stochastic  design  approach  to  various  simple  exaimple  prob¬ 
lems.  This  is  done  mainly  with  a  view  toward  illustrating  the 
improvement  in  robustness  properties  to  be  expected  and  the 
qualitative  form  of  the  control,  particularly  for  large  levels 
of  modelled  frequency  uncertainty. 

In  all  examples  considered,  the  system  retains  the  general 
form  of  (87) : 


-BK 


A  -BK-FC 
m 


(:-)-(2) 


except  that 


s 


(It 


where  represents  the  elastic  modes  while  now  comprises  both 
rigid  body  and  unstable  modes.  Matrices  B,  C,  K  and  F  are  as 
defined  previously,  and  w^  and  W2  are  white  noise  processes  with 
intensity  matrices  and  V2/  respectively. 

Uncertainty  is  modelled  only  in  the  open-loop  frequencies 
associated  with  the  elastic  modes.  It  is  assximed  that  open-loop 
frequency  deviations  are  normally  distributed  random  Vcuriables. 

In  consequence,  if  T  denotes  the  decorrelation  time  correspond- 
ing  to  the  m  elastic  mode,  then  (32)  yields: 


m 


(11' 


where  a  is  the  stamdard  deviation  of  the  m^^  elastic  mode  fre- 
quency  relative  to  its  nominal  value,  o)  .  For  illustrative  pur- 

Itl 

poses,  we  adopt  the  simple  model 


in  in 


to  reflect  a  degradation  of  structural  modelling  accuracy  with 
increasing  nominal  modal  frequency.  Thus,  in  the  following 
example  problems,  modal  frequency  uncertainty  levels  are 
uniquely  defined  by  the  relative  standard  deviation,  a,  of  the 
first  mode. 


Although  second-mean  stability  is  proved  under  the  condi¬ 
tions  stated  in  Theorem  7,  it  is  desirable  to  demonstrate  robust¬ 
ness  for  specific  designs  by  determining  stability  for  a  range  of 
parameter  variations. 

To  this  end,  we  conduct  sensitivity  studies  as  follows.  For 
each  compensator  design  (corresponding  to  particular  nominal 
values  and  a  value  of  a  in  (111) ,  we  compute  the  closed- loop  poles 
for  a  set  of  perturbed  system  models  obtained  by  replacing  the 
nominal  value  of  the  elastic  mode  dynamics  matrix,  A^,  by: 

0  l' 

A  (6)  =  block-diag  ,  _  (112) 

-2rij^a»j^(l+6)  . 

where  a  range  of  values  of  the  relative  frequency  deviation,  6, 
are  considered. 

Obviously,  since  (112)  implies  a  change  in  all  modal  fre¬ 
quencies  by  a  fixed  percentage,  robustness  for  independently 
random  frequency  deviations  cannot  be  established.  However, 
perturbations  of  the  above  form  still  provide  a  convenient  and 
practical  means  of  illustrating  relative  stability  for  a  sub¬ 
class  of  the  parameter  uncertainties  originally  postulated. 

Note  also  that  (112)  involves  variation  directly  in  the  struc¬ 
tural  mode  frequencies  as  opposed  to  parameter  variations  of  the 
form  (10) .  Thus,  although  uncertainties  in  structural  mode  fre¬ 
quencies  are  closely  allied  with  open-loop  frequency  uncertain¬ 
ties,  unconditional  stability  under  such  perturbations  may  not 
be  expected  and  (112)  provides  a  good  test  of  the  robustness 
properties  of  the  stochastic  design. 
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Further  important  features  which  we  wish  to  illustrate  in 
the  following  are  the  progressive  diagonalization  of  P  and  Q 
with  increasing  uncertainty  level  and  the  concomittant  reduction 
of  the  controller  form  to  rate  output  feedback  in  the  case  of 
colocated  actuators  and  sensors.  To  demonstrate  diagonal! zation 
of  P  and  Q,  we  introduce  the  following  measures  of  diagonal  dom¬ 
inance: 


^  =  f  -1 


Clearly,  P.,  C.  2  0,  with  equality  if  and  only  if  the  i^  column 
is  zero  except  for  the  i  element. 

4. 3  Numerical  Examples 

The  first  two  examples  discussed  here  involve  low-order 
systems  and  are  intended  to  illustrate  the  robustness  properties 
to  be  expected  under  the  stochastic  design  approach.  We  partic¬ 
ularly  consider  the  effect  of  rigid  body  and  unstable  modes  and 
non-colocation  of  actuators  and  sensors.  The  last  example  dem¬ 
onstrates  the  asymptotic  properties  of  the  control  for  large 
uncertainty  levels  for  a  fairly  high-order  system. 

A.  Two  Mass  System 

The  two  mass  systCTi  shown  in  Figure  2  provides  a  simple 
exeunple  to  Illustrate  the  relative  stability  of  the  stochastic 
design  approach.  We  assvmie  a  force  actuator  located  on  mass  1 
and  suppose  that  mass  2  is  subjected  to  a  white  disturbance 
force,  w,  of  unit  intensity.  Also,  for  simplicity,  m, ,  m^  and  k 


are  set  equal  to  unity.  Ignoring  dancing  for  the  moment,  the 
dynamic  equations  are: 


«  -Xj^  +  X2  +  u 
X2  »  Xj^  -  X2  +  w 


(114) 


It  is  desired  to  suppress  the  mean-square  displacement  of  mass  2. 
Thus,  we  choose: 


s 


E [X2  +  u^l 


(115) 


Finally,  two  cases  of  sensor  location  are  considered: 


non- colocation:  y  =  X2  +  W2 

(116) 

colocation:  y  =  *1  ^2 

where  W2  is  observation  noise  of  unit  intensity.  Note  that  in 
both  cases,  displacement  sensing  is  assumed. 

Recasting  this  probl^  in  the  modal  coordinate  basis  and 
assuming  0.5  percent  modal  damping,  the  various  matrices  appear¬ 
ing  in  (108)  and  characterizing  the  control  formulation  (88) -(90) 
become: 


"0  10  0  ■ 

O 

o 

• 

o 

0  0  0  0 

,  B  = 

0.50 

A  » 

0  0  0  1 

-0.17 

_0  0  -2  -0.014_ 

0.17_ 

Non-colocation : 

C  = 

[1/ 

o 

1 

1 

H 

Colocation: 

C  = 

[1, 

0,  1,  1] 

"  1  0 

-1 

-1~ 

0  0 

0 

0 

= 

=  1 

X 

-1  0 

1 

1 

z 

_-l  0 

1 

l_ 

O 

O 

• 

O 

0.00 

0.  00 

0.  00~ 

0.00 

0.25 

0.08 

-0.08 

,  V.  = 

X 

0.00 

0.08 

0.03 

-0.  03 

0.00 

1 

O 

• 

o 

00 

-0.03 

0.03 

a. 


b. 


1 


(117. a, b) 


(118) 


(119. a, b) 


(120. a, b) 


These  relations  define  the  nominal  system  model.  Uncertainty 
in  the  single  open- loop  elastic  mode  frequency  present  in  this 
problem  is  modelled  as  in  (110)  and  (111)  with  a  denoting  the 
standard  deviation  relative  to  the  nominal  frequency  u  =  /?  . 

Under  the  above  assumptions,  we  now  discuss  the  non- 
colocated  and  colocated  sensor  cases  separately  as  follows: 
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A. 1  Sensor  at  Mass  2  (Non-Colocation) 


Using  (103)- (106)  we  first  generated  three  designs  corre¬ 
sponding  to  uncertainty  levels  a  »  o.O,  0.5  and  1.0.  For  each 
such  design,  closed  loop  poles  were  computed  for  modal  frequency 
perturbations  of  the  form  (112) ,  with  5  varying  over  the  range 
10.90.  Note  that  since  there  is  only  one  elastic  mode,  such 
root-loci  fully  characterize  robustness  under  modal  frequency 
uncertainties . 

The  root- loci  are  shown  in  Figures  3,  4  and  5,  where  pole 
locations  for  the  nominal  system  are  indicated  by  solid  dots. 

The  deterministic  plant,  LQG  design  (a  =  0)  shown  in  Figure  3  is 
unstable  for  frequency  variations  6  =  -0.22  to  -0.90  and 
5  s  0.10  to  0.50.  It  is  seen  that  the  nominal  system  poles  at 
-0.2  ±1.4i,  corresponding  to  compensator  poles,  tend  to  push  the 
closed- loop  elastic  mode  into  the  right  half  plane.  In  contrast, 
the  stochastic  design  with  a  =  0.5  (Figure  4)  places  these 
offending  poles  at  -1.1  11. 5i  (for  the  nominal  system),  thereby 
minimizing  their  effect  on  the  elastic  mode.  This  stabilizes  the 
design  for  positive  frequency  variations  but  the  system  is  still 
unstable  for  variations  with  6  =  -0.25  to  -0.9.  When  the  model¬ 
led  uncertainty  level  is  increased  further  (a  -  1.0,  Figure  5), 
an  increased  stability  margin  is  achieved  for  positive  frequency 
variations  but  the  unstable  region  for  negative  variations 
remains  virtually  unchanged. 

A. 2  Sensor  at  Mass  1  (Colocation) 

Here,  we  again  consider  three  designs  corresponding  to 
a  8  0,  a  8  0.5  and  a  =  l.o,  and  plot  the  root- loci  for  varia¬ 
tions  of  6  from  -0.9  to  ■•■0.9.  The  results  are  shown  in  Figures 
6,  7  and  8.  In  this  case,  the  deterministic  design  (Figure  6) 
is  unstable  for  variations  6  »  -0.15  to  -0.60  and  for  6  ^  0.30. 


REAL 

TWO-MASS  SYSTEM  SENSOR  AT  MASS  2 


Fig.  3  Closed-loop  poles  for  two  mass  system, 
non-colocated  case*  a  =  0 


REAL 

TWO-MASS  SYSTEM  SENSOR  AT  MASS  2 


Pig,  4  Closed- loop  poles  for  two  mass  system, 
non-colocated  case,  o  =  0. 5 
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REAL 

TWO-MASS  SYSTEM  SENSOR  AT  MASS  2 


Fig.  5  Closed- loop  poles  for  two  mass  system, 
non-colocated  case,  a  =1.0 
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REAL 

TWO-MASS  SYSTEM  SENSOR  AT  MASS  1 


Fig.  8  Closed-loop  poles  for  two  mass  system, 
colocated  case,  o  =1.0 
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On  the  other  hand,  the  stochastic  design  for  a  =  0.5  (Figure  7) 
is  unstable  only  for  6  >  0.70.  This  stability  region  is  still 
further  extended  to  6  <  0.90  by  increasing  a  to  1.0  (Figure  8). 
Thus,  it  is  seen  that  with  sensor/actuator  colocation,  the 
stochastic  design  can  extend  the  stability  region  for  both  posi¬ 
tive  and  negative  frequency  variations.  Overall  stability  char¬ 
acteristics  may  be  shown  as  in  Figure  9,  where  we  have  combined 
the  previous  results  with  design  results  corresponding  to  addi¬ 
tional  values  of  a.  This  figure  shows  the  stability  region 
(unshaded  area)  in  the  6-a  plane.  It  is  clear  that  by  increas¬ 
ing  the  level  of  modelled  uncertainty  in  the  stochastic  design 
we  may  progressively  enlarge  the  stability  region. 

Summary;  Note  that  despite  the  presence  of  a  rigid  body  mode 
(the  effects  of  which  were  not  explicitly  treated  in  the  preced¬ 
ing  theoretical  developments) ,  no  difficulties  were  experienced 
with  the  convergence  of  (103) -(106).  Although  increase  of  the 
modelled  uncertainty  level  in  the  non-colocated  sensor  case  does 
progressively  Increase  the  overall  region  of  stability,  the 
stability  boundary  for  negative  frequency  deviations  is  not 
appreciably  affected  for  large  uncertainty  levels.  This  reflects 
an  inherent  limitation  in  non-colocated  systems  in  the  presence 
of  rigid  body  modes.  On  the  other  hand,  in  the  colocated  sensor 
case,  as  Figure  9  demonstrates,  the  region  of  stable  frequency 
deviations  increases  approximately  in  proportion  to  the  modelled 
uncertainty  level.  Thus,  it  appears  that  by  modelling  the  fre¬ 
quency  uncertainty  with  sufficient  conservatism,  an  arbitrarily 
large  stability  margin  may  be  secured. 


B.  Inverted  Pendulvim 

To  illustrate  the  effect  of  unstable  modes,  closed- loop  con¬ 
trol  of  an  inverted  pendulum,  shown  in  Figure  10,  is  investigated. 
This  system  was  previously  considered  by  Martin  [23]  and  consists 
of  a  flexible  beam  column  pinned  to  a  motorized  cart.  The  beam 
acts  as  an  inverted  pendulum  which  can  be  stabilized  by  acceler¬ 
ating  the  cart  horizontally  with  input  force  u.  Measurements  of 
the  cart  position  and  the  slope  of  the  beam  at  its  base  are 
assumed  available. 

The  system  model,  including  all  relevant  nominal  system 
parameters  are  those  given  by  Martin  [23]  and  need  not  be 
repeated  here.  It  suffices  to  note,  however,  that  the  model 
encompasses  one  rigid  body  mode,  one  unstable  mode  and  two 
elastic  modes.  The  unstable  mode  frequency  is  ±3.22  rad/sec 
while  the  first  two  bending  mode  frequencies  are  3.055  and 
22.05  rad/sec.  The  inherent  damping  factor  for  the  bending 
modes  is  assvimed  to  be  0.  005.  Finally,  the  state  weighting 
matrix  is  diagonal: 


=  diag  [0.001,  0,  0.4,  0.1,  0,  0.1,  0] 


(121) 


Uncertainties  only  in  the  bending  mode  frequencies  are  con¬ 
sidered  here  and  these  are  modelled  in  accordance  with  (110)  and 
(111) . 

With  these  assumptions,  three  designs  were  computed  corre¬ 
sponding  to  a  =  0,  0.2  and  0.5,  where  a  is  the  first  elastic 
mode  standard  deviation.  Root- loci  for  elastic  mode  frequency 
deviations  of  the  form  (112)  with  6  in  the  range  ±0.30  are 
shown  in  Figures  11  through  13.  Notice  that,  in  the  deterministic 


77 


Y////////77//7//Z7/7/. 

Fig.  10  Inverted  pendulum  system 
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IMAGINARY  PART 


design  (a  s  0,  Figure  11)/  both  bending  inodes  become  unstable 
for  |5|  >  0.03.  Increase  of  the  modelled  uncertainty  to  a  =  0.2 
(Figure  12)  pushes  the  compensator  poles  deeper  into  the  left 
half  plane  and  dramatically  increases  the  stability  region  to 
6  et-0.09,  0.29].  A  still  larger  value  of  a  (a  =  0.5,  Figure  13) 
completely  stabilizes  the  second  bending  mode  over  the  range  of 
frequency  deviations  considered.  On  the  other  hand,  while  the 
first  bending  mode  is  also  stabilized  for  all  positive  frequency 
variations,  it  becomes  unstable  for  6  <  -0.1.  In  other  words. 
Increasingly  large  levels  of  modelled  uncertainty  result  in 
relatively  minor  increases  in  the  stability  margin  for  negative 
frequency  variations. 

Summary ;  Although  the  theoretical  developments  of  earlier 
chapters  assumed  a  stable  open-loop  system,  the  computational 
scheme  of  (103) -(106)  converged  without  difficulty  in  this  case 
to  produce  stable,  robust  designs.  However,  the  above  results 
do  indicate  that  the  presence  of  an  unstable  mode  limits  the 
ability  of  the  stochastic  design  approach  to  improve  the  stabil¬ 
ity  margin  for  negative  frequency  deviations  -  particularly  for 
the  elastic  mode  in  closest  proximity  to  the  unstable  mode.  In 
part,  this  may  be  ascribed  to  the  fact  that  in  modelling  uncer¬ 
tainty  only  in  the  elastic  modal  frequencies  we  ignore  the  influ¬ 
ence  of  such  uncertainties  on  the  unstable  mode  so  that  robust¬ 
ness  improvement  is  mainly  restricted  to  those  modes  (i.e.,  the 
second  bending  mode,  in  this  case)  which  are  accurately  repre¬ 
sented  in  the  stochastic  design  model.  Nevertheless,  the 
stochastic  design  approach  yields  an  enormous  enlargement  of  the 
stability  region  over  that  obtained  by  the  conventional  LQG 
design. 


Here  we  consider  a  simply-supported  Bernoulli-Euler  beam 
with  normalized  span-wise  coordinate  as  shown  in  Figure  14.  A 
force  actuator  and  colocated  rate  sensor  are  assumed  at 
5  =  C  =  2/43.  The  nominal  system  dynamics  matrix  in  the  modal 
coordinate  basis  is  block-diagonal  with  diagonal  blocks  of  the 
form  (3) .  With  appropriate  non-dimensionalization  of  the  equa¬ 
tions  of  motion,  we  may  write: 


0) 


k 


k  =  1, . . .  ,n 


(122) 


Furthermore,  we  set: 


Tij^  =  0.  005  ,  V  k 


(123) 


and  B  and  C  assume  the  form: 


B  = 

®kl  ’  I  (124) 

k  =  1, . . . , 2n 


For  this  exeunple,  an  "energy"  state  weighting  is  chosen; 


Fig.  14  Simply- supported  beeun  with  colocated  force 
actuator  and  rate  sensor 
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0 


(125) 


R 


1 


block-diag 
k=l, . . . ,n 


xj 


and  we  suppose  the  beam  to  be  excited  by  a  spatially  and  tempo¬ 
rally  white  disturbance  force  with  unit  intensity,  i.e.,; 


V 


1 


0  0 

block-diag 

k=l, . . . ,n  0  1 


(126) 


Since  only  one  actuator  and  sensor  are  used,  R2  and  V2  are 
positive  scalars  which  we  denote  by  and  p^,  respectively. 

As  before,  open-loop  frequency  uncertainties  are  modelled 
according  to  (110)  and  (111)  so  that  uncertainty  levels  are  com¬ 
pletely  specified  by  the  standard  deviation,  a,  of  the  first  mode. 
Clearly,  (111)  also  entails  a  fairly  rapid  increase  of  uncer¬ 
tainty  level  with  increasing  modal  order. 

Under  the  above  conditions,  stochastic  designs  were  com¬ 
puted  for  a  range  of  values  of  a,  p^  and  p^.  For  each  such 
design,  root- loci  were  determined  for  system  perturbations  of 
the  form  (112) .  It  suffices  to  note  that  both  stochastic  and 
deterministic  designs  were  stable  for  all  frequency  variations 
considered.  Thus,  improvement  in  robustness  is  not  an  important 
issue  in  this  example  problem. 

What  this  exeunple  shows  most  clearly  is  the  distinctive 
form  of  the  control  provided  by  the  stochastic  approach. 

Figures  15  through  20  contrast  the  deterministic  design  (a  =  0) 
with  a  stochastic  design  (a  =  0.2),  both  with  p^  0.1  and 
»  0.1  and  15  modes  retained  in  the  design  model. 


I 

Figures  15  and  16  show  the  diagonal  dominance  indicators 

!  defined  by  (113).  Recall  that  P.,  Q,  =  0  implies  that  only  the 

i  tH  ^  ^ 

I  diagonal  element  of  the  k  row  of  P  and  Q,  respectively,  are 

non- zero.  Consequently,  the  results  for  the  stochastic  design 
show  that  except  for  the  sub-blocks  corresponding  to  the  first 
~8  modes,  P  and  Q  are  nearly  diagonal.  In  contrast,  the  indi¬ 
cators  Pj^  and  for  the  deteanninistic  design  remain  signifi¬ 
cantly  above  zero  for  all  modes  considered. 

Diagonalization  of  P  and  Q  in  the  stochastic  design  results 
in  suppression  of  displacement  loops  as  can  be  seen  from  the 
regulator  and  observer  gains  displayed  in  Figures  17-20.  While 
regulator  position  gains  increase  with  modal  order  in  the  deter¬ 
ministic  design  (Figure  17) ,  the  stochastic  position  gains  are 
negligible  beyond  the  8th  mode.  At  the  same  time,  velocity 
gains  (Figure  18)  are  nearly  the  same  in  both  designs.  The  same 
general  tendencies  can  be  seen  for  the  observer  gains  (Figures 
19  and  20) .  Thus,  displacement  loops  are  almost  completely 
suppressed  for  the  higher-order,  relatively  more  uncertain  modes 
and  the  solution  for  the  corresponding  portions  of  P  and  Q 
reduces  to  the  asymptotic  forms  given  by  Theorem  8. 

It  should  also  be  noted  from  Figures  17-20  that  gains  for 
the  first  three  modes  (corresponding  to  more  than  20%  closed-loop 
deunping  for  these  modes)  are  nearly  the  same  for  both  determinis¬ 
tic  and  stochastic  designs.  This  indicates  a  region  of  "quasi- 
deterministic"  control  for  low  order  modes. 

Summary  t  This  case  conforms  to  the  conditions  assumed  in  the 
theorems  of  Chapter  3  and  it  is  no  surprise  that  robust  stability 
is  readily  achieved.  Moreover,  it  is  evident  that  the  stochasti¬ 
cally  designed  control  for  modes  beyond  the  8th  closely  approxi¬ 
mates  the  asymptotic  form  given  by  Theorem  8  and  Lemma  8.  In 
fact,  since  the  asymptotic  control  is  rapidly  approached  with 
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Fig.  15  Diagonal  dominance  indicator  for  regulator 
cost  matrix 
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Fig*  16  Diagonal  Dominance  indicator  for  observer 
cost  matrix 


88 


flD-fll31  276  MEflN-SQUflRE  OPTIMAL  FULL-ORDER  COMPENSATION  OF  2/2 

STRUCTURAL  SVSTEMS  WITH  UN.  .  CU)  MASSACHUSETTS  INST  OF 
TECH  LEXINGTON  LINCOLN  LAB  D  C  HVLAND  01  JUN  83  TR-626 
UNCLASSIFIED  ESD-TR-83-013  F19628-80-C-0002  F/G  12/1  NL 


Fig.  18  Velocity  gains  for  stochastic  and  deterministic 
regulator 


ABSOLUTE  VALUES  OF  POSITION  GAINS  FOR  OBSERVER 


BEAN  15  NODE  DESIGN 


I  Z  3  H  5  6  7  6  9  10  1  112  13  14  15 

MODE  NUMBER 


LEGEND 


S06PGN 

OOBPGN 


0.2 

0.0 


Fig*  19  Position  gains  for  stochastic  and  deterministic 
observer 
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increasing  modal  order,  (93)  may  be  used  to  determine  the  gains 
for  as  many  modes  beyond  the  15th  as  desired  with  no  significant 
computational  burden.  Because  of  the  properties  noted  above  for 
the  control  of  low-order  modes,  we  may  say  in  summary  that  the 
stochastic  design  approach  automatically  produces  a  high  author¬ 
ity,  essentially  deterministic  design  for  low-order,  well-known 
modes  and  a  low  authority,  rate  output  feedback  control  for  high- 
order,  very  uncertain  modes.  These  two  regimes  are  seen  to  exist 
as  limiting  qualitative  features  of  a  unified,  globally  stable 
design. 


5. 


CONCLUSION 


In  the  foregoing  developments,  the  minimvim  information 
approach  to  the  modelling  of  uncertain  structural  systems  has 
been  applied  to  the  mean-square  optimal  design  of  full-order 
dynamic  compensation.  It  was  seen  that  judicious  choice  of  the 
coBQ>ensator  dyn2unics  matrix  permits  the  optimality  conditions 
to  be  reduced  to  relatively  simple  forms.  Under  fairly  mild 
restrictions,  these  optimality  conditions  were  shown  to  possess 
unique  solutions  for  which  closed-loop  stochastic  stability  is 
quaranteed.  Furthermore,  the  control  reduces  to  a  simple, 
inherently  robust  asymptotic  form  for  sufficiently  high  levels 
of  modelled  uncertainty. 

A  simple,  straightforward  computational  scheme  was  devised 
for  numerical  solution  of  the  optimality  conditions.  The  first 
two  numerical  examples  considered  in  the  last  chapter  not  only 
show  this  algorithm  capable  of  handling  rigid  body  and  open- loop 
unstable  modes  but  also  suggest  that  the  conditions  assumed  in 
the  existence  and  uniqueness  result  of  Chapter  3  can  be  consid¬ 
erably  relaxed.  Moreover,  these  example  problems  show  that  the 
stochastic  design  approach  can  enormously  increase  stability 
margins  over  what  can  be  achieved  by  a  conventional,  LQG  design. 

Furthermore,  the  last  numerical  example  discussed  in 
Section  4.3  suggests  that  when  reciprocal  decorrelation  times 
increase  monotonically  with  modal  order,  closed- loop  modes  may 
be  divided  into  two  qualitative  regimes:  the  "coherent"  and 
"incoherent"  systems  (to  use  the  terms  introduced  in  Ref.  [7]). 
Incoherent  modes  are  associated  with  great  a  priori  uncertainty 
and  are  mutually  uncorrelated  and  uncorrelated  with  the  coherent 
system  composed  of  low-order,  relatively  well-known  modes.  Since 
the  mean-square  optimal  control  for  the  Incoherent  system  is 
known  in  advance  (by  virtue  of  Theorem  8  and  expressions  (93)), 


these  qualitative  features  permit  the  solution  of  (89)  and  (90) 
for  very  high-order  systems  by  combination  of  (93)  with  the  solu¬ 
tion  of  the  reduced-order  forms  of  (89)  and  (90)  for  the  coher¬ 
ent  syston. 

Such  a  scheme  would  reduce  the  computational  burden  to  that 
associated  with  the  relatively  few  well-known  modes.  Of  course, 
this  would  obviate  the  difficulties  of  design  computation  but  not 
of  implementation.  Although  the  asymptotic  form  of  control  for 
the  incoherent  system  approximates  a  rate  output  feedback  law, 
thereby  aiding  practical  implementation,  a  completely  satisfac¬ 
tory  treatment  must  await  the  extension  of  the  stochastic  design 
approach  to  fixed-order  dynamic  compensation. 
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APPENDIX  1 

Proof  of  Theorem  4 

First,  we  write  (46. c)  and  (49)  in  accordance  with  the  state 
transformation  (50).  Using  (51),  and  the  notation: 


A 


(A. 1.1) 


equations  (46. c)  and  (49. a)  yield 


% 

1 


(V-6K)«  Pj  .  c"?!,  . 

.  r,  .  k«R,.  .  r(P^  .  P^^  .  pH^  .  P^) 


P^q  +  c“Pq  +  P^Bk  +  +  «)  b.  >  (A. 1.2) 

-  <"^2'^  I 


I 


K®6“p^g  +  (6ic  +  a)“  Pg  +  p“q6K  +  Pq(6K  +  a)  c. 
+  ic**R2K 


+  Q^qK  B"  +  +  Vj^ 


(Mj^-Bic)  Q^q  +  Bk  Qq  +  Q^c”  +  +  Sk)“  b. 


(A. 1.3) 


+  J{Q^}  + 


0  =  CQ^q  +  (a  +  Bk)  Qq  +  Q^g?”  +  Qg(“  +  c. 


+  l{Q^}  +  +  f 


and  (49.b,c)  become: 


-Oeq  *  °q> 


=  Rj6  1^(VV  *  '‘5q‘°“q'V 


(A.1.4) 


^  *wll  000  e^tt  11 


(A. 1.5) 


while  the  identity,  (49.d) ,  assumes  the  form: 


iq<05q-°C>  ^  ^'V°Cq>  ’  ^ 


(A. 1.6) 


Now  taking  C  ^  0  (i.e, ,  choosing  a  as  in  (53)),  we  show  that 
(54)  through  (57)  identically  satisfy  (A. 1.2)  through  (A. 1.6). 

First  note  that  with  =  0,  (A. 1.6)  becomes 


-°1q' 


0 


which  is  identically  satisfied  by  (54. b) .  With  (54),  (A. 1.4)  and 
(A.  1.5)  assume  the  foms 


[K  Pjl  (Qj  -Qq)  -  0 

R„[f  -sy  v:h  -  0 

q  q  2 

and  these  are  satisfied  by  (55).  By  virtue  of  (55.b),  it  is 

readily  checked  that  with  ?  -  0,  (A.1.3.b)  and  (A.1.3.C)  yield 

Q  Q?  in  accordauice  with  (54.b).  Thus,  Q>  and  equation 
q  5q  ^<3 

(A. 1.3)  need  not  be  considered  further.  Moreover,  with  (54. a) 

and  (55. a),  (A.1.2.b)  becomes: 


0  .  Pq 


which  is  satisfied  if  ^  vanishes. 


In  consequence,  only  (A. 1.2. a),  (A.1.2.c),  (A. 1.3. a)  and 
(A.1.3.C)  remain  for  consideration.  With  C  =  0  and  (54)  these 
equations  become  (56)  and  (57). 

Thus,  Eqs.  (54)  through  (57)  identically  satisfy  the  sta¬ 
tionary  conditions  (46. c)  and  (49). 
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APPENDIX  2 


Proof  of  Theorem  5 

Since  only  elastic  modes  cure  considered  and  n,,  >  0  for  all 
—  h  h  — 

k,  it  follows  that  (y,a  )  is  stabilizable  and  (r?,y)  is  detect- 

/V 

ed>le.  As  P  is  presumed  to  exist  as  a  positive  semi-definite 
matrix,  ( (r +  j{p})^,il)  is  also  detectable  by  Theorem  3.6  of 
Ref.  [20].  In  other  words  (62. a)  assumes  the  form  of  the 
stochastic  Riccati  equation  for  the  regulator  problem: 


0  =  vi“p  +  Py  +  J{P}  +  s  -P  a  P 

m  m 

s  =  rj^  +  j{p}  2  0 

—  j, 

(y,  a  )  stabilizable 
h  — 

(s  , y)  detectable 


(A. 2.1) 


Then  by  virtue  of  Theorem  12  of  Ref.  [7],  y  -ap  is  asymptotically 
stable.  But  -1/2  is  non-positive  and  diagonal  so  that  y^  -aP  is 
also  asymptotically  stable.  Analogous  reasoning  for  (62.b) 
establishes  the  asserted  stability  properties  of  (y  -Q  o)  and 
()]^  -Q  o)  and  completes  the  proof  of  part  A. 

Consider  system  (58)  in  the  absence  of  disturbance  and 
observation  noise: 


5  =  yC  -BKq 


(A. 2. 2) 


q  =  (Pjjj  -Bic)q  +  fY(?-q) 


with  <  and  f  given  by  (61).  To  prove  part  B  we  show  exponential 
asymptotic  stability  of  the  second  moment  response  of  this  system. 
Partitioning  the  second  moment  matrix  of  in  accordance  with 


(48),  so  that 


Q  =  E 


10 


(5®,  q”) 


.H 


Q 


Kq 


5q 


the  Lyapunov  equation  for  Q  yields: 


Q, 


“mOs  -6''  Olq  ^ 


''mOEq  °q  *  *  °Eq'“m  -‘if 


'Cq 


H 


H 


Q^q  (l'm“^Y-B'c)Qg  +  Q^gyV  +  QqCy^-fY-BK)  -  c. 


Alternately,  letting: 


5  ^ 


°E 

■-°5q 


Q 


?q 

/W 

Q_ 


=  T  Q  T 


(A. 2. 3) 


(A. 2. 4) 


(A. 2. 5) 
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with  T  defined  by  (50) ,  we  have: 


°q  =  <%-*1f>0q  *  ^ 


Assuming  P  and  P  to  be  positive  semi-definite  solutions 
(62)  and  (63)  ,  consider  the  non-negative  quantity: 


£ 


A 


0 

P  - 


Q 


Use  of  (62. a),  (63. a)  and  (A. 2. 7)  yields: 


£  =  -tr[Q^r^  t  tl2n'-"2nJ 


Q? 

?q 


"5q 


][.: 


2n 


2n 


POP] 


=  -trlQ^rj^  +  Qq  P  0  P]  <  0 


where  the  last  line  follows  from  (A. 2. 5).  Since  Q_P  a  P  is 

H  H  ^ 

simply  E[(^k  R2Kq],  we  conclude  from  (A.  2. 8)  that  for  all 

«  Q(t»0)  ^  0,  6Kq  converges  exponentially  to  zero  in  the 

mean-square  as  t  increases  without  bound.  Furthermore: 


(A. 2. 6) 


to 


(A.2.7) 


(A. 2. 8) 


(A. 2. 9) 


|(eKQ”q)j,j|  ^  ((E(6Kq  q“ic“B“))j,j,(E(5 

In  view  of  (A.  2. 9)  and  because  A^-[u  +  Qr  y  +  is 

QC  in  ^  ^  in  ^ 

asymptotically  stable  (see  proof  of  TheorCTi  8,  Ref.  [7]), 

(A. 2. 4. a)  shows  that  must  converge  exponentially  to  zero. 

Now  consider  (A.2.6.b,c).  The  term  J{Q^}  in  (A.2.6.C)  con¬ 
verges  exponentially  to  zero  while  is  asymptotically 

stable.  Thus,  the  solution,  of  (A.2.6.C)  converges  exponen¬ 
tially.  Since  both  (y^-BK)  and  are  asymptotically  stable, 

similar  reasoning  applied  to  (A.2.6.b)  shows  that  also  con¬ 

verges  exponentially.  Consequently,  all  elements  of  Q  converge 
exponentially  to  zero  as  t  increases  without  bound  and  system 
(A. 2. 2)  is  second  mean  stable.  Almost  sure  exponential  stability 
is  then  directly  implied. 

Finally,  since  (61)  through  (63)  imply  the  stationary  con¬ 
dition  and  the  control  given  by  (61)  is  admissible,  the  conclu¬ 
sion  of  part  C  follows  by  Lemma  2. 


APPENDIX  3 


As  a  preliminary  step  to  the  proof  of  Theorem  6  we  show: 
Lemma  A. 3. 1 

Under  the  conditions  of  Theorem  6  suppose  that  P^,  >  0. 

Then: 


controllable 

irl.vj  , 

reconstructible 

((PiCrt’i)’*,  Iv-Qii) 

reconstructible 

(y  -aP.,  (Q.aQ.)^) 
in  X  XI 

controllable 

Proof 

_  j, 

(li,a  )  is  controllable  if  and  only  if  (see  Ref.  [20],  p.  45) 

rank  [p  -XI,  a  ]  =  2n 

for  all  X  defined  on  the  spectrum  of  y.  Since  y  and  I  are  both 
diagonal  and  the  modal  frequencies,  u)^,  k  »  l,...,n  are  assumed 
distinct,  this  implies: 

— 

rank  Iy„  -XI,  a^]  =  2n 

nl 
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—  —  k 

for  all  X  on  the  spectrum  of  y_.  Therefore  (y  ,a^)  is  control- 

®  —  if  “J*  —  —k  — 

lable.  The  stated  conclusions  for  (y«.»v7)  ,  (r7,y_)  and  (^,y_,) 

ml  1  m  m 

can  be  shown  in  a  similar  manner. 

Now  consider  the  matrix  [y  -Xl>P.S]  with  X  on  the  spectrum 

T|  ^ 

of  vr^.  Writing  this  in  the  modal  coordinate  basis  (in  which  the 
m  ’  /  «  -1  \ 


/  0  1  \ 

nominal  syst&n  map,  y,  has  the  matrix  block-diag  I  _2  _  1 ) 

and  noting  that  (y**,6)  is  ~^'^k‘**k 

controllable,  it  is  seen  that  the  rank  of  [y  -XI, P.8]  is  less 
than  2n  only  if  6  is  orthogonal  to  more  than  one  row  of  P^.  But 
since  >  0  this  is  impossible.  Thus. 


rank  [y^  -XI,  p^  6] 


for  all  X  on  the  spectrum  of  y^.  This  implies: 


2n-l  „  . 

2  ’'p  =  c  " 

k=0  “ 


where  p  denotes  the  image  of  P.B.  In  consequence,  denoting  the 

H  ^ 

image  of  y  by  (?: 


i.e.,  (y^  -oQj  B)  is  control l2a>le.  Since  >  0  this  means 
_  m  1  ji  2 

(Pjij  -oQ^,  (P^oPj^P  )  is  also  controllable.  Therefore 

( (P.  op. ) o)  is  reconstructible.  In  exactly  analogous 
XX  in  X  _  L 

fashion,  (y  -oP, , (Q4 oQ. ) ^)  may  be  shown  to  be  controllable, 
in  X  XX 

The  following  helps  to  show  that  all  members  of  the 
sequences  defined  by  (67)  and  (68)  can  be  made  positive  definite: 

Lemma  A. 3 . 2 

Under  the  conditions  of  Theorem  6,  P  >0  and  Q  >0 

o  o 

together  with  (67)  and  (68)  imply: 


>  0  ;  i  i  0 

(y  -oP.)  &  (y  -Q.o)  stable  ,  i  >  0 
in  X  m  X 

^  A 

Pi»  Qi  >  0  ;  i  2  1 

Proof 

^  A 

Suppose  that  P.  >  0  and  >  0.  Then  by  Theorem  3.6  of 
Ref.  [20)  (  +  i’{Pj^)) **,y)  and  (jT,  (Vj^  +  i{Qj^})**)  are  recon¬ 

structible  and  controllable,  respectively.  By  virtue  of  Lemma  6 


and  (67) ,  p^^  auid  are  positive  definite  and  (ym-^Pi)  and 
(y^-Qj^o)  are  asymptotically  stable. 


by  virtue  of  L^nma  4. A  and  (68),  and  are  positive  defi¬ 


nite.  The  stated  conclusion  then  follows  by  induction  on  i.  □ 


As  a  final  preliminary  step,  we  specifically  establish 
quantities  which  will  serve  as  upper  and  lower  bounds  for  the 
sequences  defined  by  (67)  and  (68). 

Lemma  A. 3. 3 

Positive  definite  p'^,  p'*,  are  uniquely  defined  by: 


u—  u  u 

0  *  p  u  +  ir^P  +  P  op 
m  xn 


a. 


0  =  y^p'^  +  p'^y.  +  r.  +  Kp'^}  +  I{V^}  -P^aP'^  b. 

m  m  1 


(A. 3.1) 


0  -  Q“*5  +  +  Q"^Q“ 

0  -  ivo"  +  Q"iS  +  »!  +  «Q“>  + 

xn  xn  1 


a. 


b. 


(A.3.2) 


Likewise,  the  equations: 

0  =  P^(y  -Q^o)  +  (y  -Q“a)“p^  +  P^aP^ 
xn  in 

» - *  -^1  * 

0  -  Q^(y„-ap'^)“  +  (y^-ap'^)Q*’  +  Q*'oQ^ 
0  “  y«Q*’  +  +  V,  +  I{Q^}  -Q^oQ^ 


a. 


b. 


I  (A. 3. 3) 


a. 


b. 


(A.3.4) 


'^1  "^0 

uniquely  define  positive  definite  P,  P  #  Q  and  Q  with: 


P^  ^  p'^  f  P^  ^ 


^u  u 

Q  ^  Q  f  Q  ^ 


I  (A.  3. 5) 


Moreover,  ‘  stable. 


Proof 


First  consider  (A. 3.1).  From  (A. 3.1. a)  and  the  assumption 


Xj^(y)  <  0: 


liP^}  <{P%P"} 


Then,  by  Lemma  6  and  (A. 3.1.b)  it  follows  that  P^  ^  P*  where: 


“  ***’^m  "*■  ^1  ■*■  ^P*oP*>  -P*crP*  (A. 3. 6) 


Hence  the  diagonal  portion  of  P*  satisfies; 


0  -  y"{P*}  +  {P*}y  +  {r^^} 
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As  X  (y)  <  0,  ^  0,  this  uniquely  defines  0  s  {P*}  <  ®. 

R  X  ^ 

Since  P*  is  hermitian,  it  follows  that  P*  and  hence  P  is  bounded. 

Likewise,  from  (A. 3.1. a),  p“  is  bounded  and  positive  definite. 

With  the  aid  of  Lemmas  4  and  6,  the  sequences  pV  and  pV 
(i  2  0)  defined  by 


0 

0 


y**P?  +  pVy  +  r,  +  I{pV  +  pV}  -P?aPY 
^m  1  i^m  1  1  1  11 


(A.3.7) 


may  be  shown  to  be  positive  definite,  monotone  nondecreasing  and, 
by  the  previous  argument,  bounded  from  above.  Use  of  Lemma  3 
then  implies  that  (A. 3.1)  possesses  positive  definite  solutions 
for  p'^  and  p'^.  To  show  uniqueness,  let  (P^,  P^)  and  (p^,  P^) 
two  sets  of  solutions.  Manipulation  of  (A. 3.1)  produces: 


0  =  y”z  +  Z  y„  +  J{Z} 

m  iR 

where  Z  denotes  (P^  +  P^  -P^  “^2^ •  It  is  easily  seen  that 

A- Iy?Z  +  ZiT  +  I{Z}]  is  asymptotically  stable,  whence  Z  =  0 
Z  m  m  u  ''u 

uniquely  by  Lonma  5.  Thus  P  -f  P  in  unique,  and  (A.3.1.b) 

assumes  the  form: 


?!  -  +  p“}  a  0 

where  is  uniquely  defined  and  (r^#  is  reconstructible. 

In  consequence,  the  above  Riccati  equation  uniquely  defines 
I  as  a  positive  definite  matrix. 

Thus  (A. 3.1)  possesses  unique  positive  definite  solutions 

M 

for  P  and  P  .  Moreover,  application  of  Lemma  6  shows  that 
(iTjj  -aP  )  is  asymptotically  stable.  The  stated  results  for 
(A. 3.2)  follow  analogously. 

Finally  consider  (A. 3. 3).  Application  of  Lemma  6  to 
(A.3.3.b)  and  Lemmas  4  and  A. 3.1  to  (A. 3. 3. a)  suffice  to  show 
the  existence  and  uniqueness  of  positive  definite  solutions  and 
the  stability  of  (y^  -aP^) .  The  results  for  (A. 3. 4)  follow 
similarly. 

>s 

Note  that  since  >  0,  comparison  of  (A.3.1.b)  and 

D  M 

(A. 3.3.b)  yields  P  ^  P  by  Lemma  6.  In  consequence,  use  of 
Lemma  4  on  (A. 3.1. a)  and  (A. 3. 3. a)  suffices  to  show  P^  ^  P^. 
Analogous  proof  of  the  remaining  properties,  (A. 3. 5),  is 
straightforward.  □ 

Mow  we  are  in  a  position  to  prove  Theorem  6. 

Proof  of  Theorem  6 
Suppose  that 


i  P 


0  < 


(A. 3. 8) 
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From  the  previous  Lemma  P  >  0  so  that  >  0  by  Lemma  A.  3.2. 
Hence,  Lanma  6  and  comparison  of  (A.3.3.b)  and  (67. a)  yields: 


P. 

1 


>  P 


i 


Similarly: 


Qi  > 


Also,  since  <  Q  ,  it  follows  that 


Q-  s 
1 


nU 


from  Lemma  6  and  comparison  of  (A. 3.4.b) 

P^  s  p". 

Manipulation  of  (68. a)  and  (A. 3. 3. a) 


and  (67.b) . 

gives: 


Similarly, 


0  =  (Pi^.1  -P^)  (iTj^-Q^a)  +  (P.^i  -P*-) 


i+1 


/V  /\  i 


(A. 3. 9) 


K  0  ^  /K  0 

+  (P^  -P*)  a  (P^  -P^)  +  S 


where 
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i 


s  ^  Pi+i<Q“  -Qi)  o  +  a  (q'^ 


(A.  3. 10) 


+  P*'aP^  +  P^OP^ 


Now  -Q.  s  0  and  a  ^  0  and,  since  P..,  >  0  by  Lemma  A.  3.3, 

^  —  —  u 

the  positive  semi-definiteness  of  -Qj^)  +  ^(Q 

follows  by  use  of  Lemma  4.  Similarly,  (P^aP^  +  p^ap^)  ^  0  so 

that  S  ^  0.  Also,  because  P^  -P*"  >  0,  (( (P^^  -P^)a(P^  -P*")  )^, 

y_  -Q^o)  is  reconstructible  by  Lemma  A. 3.1.  Then,  since 

(y^  -Q  o)  is  asymptotically  stable.  Lemma  4. A  gives 


^+1  >  P 


Repetition  of  the  same  sort  of  argument  yields 

A  /V0  /V 

Oi+1  ^  ^  ^  and  ^  Q  •  In  summary,  we  have  shown 

that  (A. 3.8)  implies: 


P^  <  P.  s  <  Q.  ^ 


P  <  P.^i  i  P"  ,  Q*-  <  Q,^,  s  o'" 


Induction  on  i  completes  the  proof. 
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